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ABSTRACT 



The purpose of this thesis is the investigation of 
mathematical methods for treating collisions of light nuclei. 

Our understanding of nuclear forces is derived largely 
from the study of simple nuclear systems. Bound nuclei 
consisting of two to four particles have been theoretically 
treated with considerable accuracy and yielded many valuable 
results. Another important source of information are nuclear 
collisions. However, up to the present, only two particle 
collisions have been theoretically treated in a systematic 
and reliable manner. 

This thesis is concerned with collisions involving any 
number of nucleons. The problem is to calculate the scattering 
properties of nuclei from assumed nuclear interactions. These 
theoretical results can then be checked against experimental 
data, thus providing a test for the validity of the interactions 
which were postulated. As the wave function of a nuclear 
n - body system extends over 3n dimensions and is highly 
irregular, on account of the great strength of nuclear forces, 
one cannot hope to determine it with great exactitude for 
systems involving more than two particles. For this reason 
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we have developed variational methods which yield the 
scattering properties of nuclei with considerable accuracy 
without requiring an exact knowledge of the wave function* 

The wave function is assumed to satisfy a Schroedinger 
equation in which nuclear forces are described by two-particle, 
velocity-independent potentials* When the colliding nuclei 
are far apart, the form of the wave function is known: It 

consists of a product of functions describing the internal 
states of the nuclei, and a free particle function, F(r ), of 
the relative position vector r* of their centers of gravity. 

To predict the experimental scattering properties of the 
nuclei the detailed characteristics of F(lr ) must be 
determined from the assumed interactions. 

If one considers collisions in which the nuclei have 
a definite angular momentum of relative motion , Z , the 
function F( t* ) ha3 the asymptotic form 

£ £ f ) — — > ^ S Via il A rj ^ Xz /v ^ S ^ 

where k Is the de Broglie v/ave number of relative motion. 

As. \ y \ , is the angular function and It the phase 
shift corresponding to the angular momentum • The phase 
shift determines the scattering cross section. 

Similarly, the function F(P ) corresponding to a 
plane wave of one nucleus incident on another has the 
asymptotic form 
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Here ( k ( | >- ) k I •=• k, k, is in the direction of incidence 

and f(k | , k) is the scattering amplitude in the direction 
of k. 

We seek to determine the phase shifts and scattering 
amplitudes, given the interaction between the nucleons. In 
collisions of two nucleons the wave equation separates and 
its solution presents no essential difficulty. But when the 
colliding nuclei are complex, the Schroedinger equation can 
no longer be separated. Hence an exact solution is not 
possible and one must take recourse to approximation methods. 

We have followed tv/o lines of attack, both of which 
have been successful in bound state problems: The Rayleigh- 

Ritz variational method of undetermined coefficients and the 
variation-iteration method. 

Raylelgh-Rltz Method * 

The principle of one of the procedures adopted may be 
illustrated by the one dimensional Schroedinger equation, 

| -t 4 - \/(*)j ^(.x) - o 

which is to be solved subject to the boundary condition that 
u(0) -s. 0. The interaction potential V(x) may be supposed to 
vanish for x ^ a, so that 

5 S* i/H A *V A X, } DC > CL 



We now form the expression 
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and consider its first variation, restricting our trial 
functions, u^, , by the conditions u^fO) =• 0 and 

(*) * — ^ "?C "V ^ Csyo X x — J) q T> 

After an integration by parts, we find 

5 J rs. - 4 ^ ^ ^ : ( x ■+ JL • ^ o 

Since for the correct solution u, I ^ 0, it follows that 
tan determined from the equation 

Jk 'tow n ■= X 

is stationary, relative to any admissible variation of u. 

When this equation has been re-written in a homogeneous 
form, one can use a trial function of the type 

i w i ^ ^2. ^ Z ^ ^ h ^u. 

on the right hand side and determine the c's from the 
stationary property of tan \^ . As in the case of bound 
state problems this leads to a de terminantal compatibility 
equation. In contrast to secular equations this is a linear 
equation for the determination of tan vj , regardless of the 
number of parameters. 

This and similar methods are applied to the following 



cases : 
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A. Two Particle Collisions. 

1. Phase Shifts. Stationary expressions are developed 
for the phase shifts in the case of central forces and for 
the proper phase shifts of the scattering matrix in thB case 
of tensor forces. Some numerical calculations show that an 
accuracy of a fraction of one per cent is attainable without 
undue labour. 

2. Scattering Amplitudes. The variational method 
developed permits -an application of the Rayleigh-Ritz method; 
it also leads naturally to the Born approximations and 
Schwinger's variation-iteration formulation ( J.S. Schwinger, 
unpublished lectures ). 

3. Level Width. A stationary expression for the width 
of a resonance level is derived. 

B. Collisions of Composite Nuclei. 

1. Neutron Deuteron Scattering. This case is treated 
in some detail, as a typical example of composite collisions. 

A preliminary calculation is carried out for the scattering 
cross section at zero energy and reasonable agreement 

( within about ten per cent ) is obtained with the experimental 
value of 3.2 barns. 

2. General Composite Collisions. A few features of 
the Rayleigh-Ritz method in rather general collisions, which 
may involve excitation or transmutation, are discussed. 
Variational principles for the elements of the scattering 
matrix are derived. 
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Variation Iteration Method . 

Instead of employing trial functions with undetermined 
coefficients in the stationary expressions for the phase shifts 
and scattering amplitudes one can also use iterated functions. 

Thus in the one dimensional problem the n-th iterate 
is defined in terms of the (n-l)-st by the equation 



Professor Schwinger ( J. S. Schwinger, unpublished lectures ) 
has applied this method very successfully In collisions of 
two particles. Since the iteration improves the wave function 
very rapidly this method converges extremely well. Schwinger’s 
method is generalized from two particle collisions to 
composite collisions. The persistence of the internal 
interactions in the asymptotic part of the configuration 
space causes considerable difficulty In the construction of 
a suitable Green’s function. Of several methods which are 
investigated only one ( involving a Green’s function in three 
dimensions ), is found at all feasible. However even this 
procedure is so complicated, that unless a more suitable 
version of the variation-iteration method can be found, the 
Rayleigh-Ritz method appears to be more useful for the treatment 
of composite nuclear collisions. 
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INTRODUCTION. 



A great deal of information about the nature of nuclear 
forces ha 3 been derived from a comparison between experimental 
and theoretical studies of simple nuclear systems. Among the 

most important data so compared are the binding energy and the 

1-1 

magnetic dipole and quadrupole moments of the deuteron ; the 

^ 3 4 

binding energies of H , He and He ; and the scattering 
cross sections of neutrons and protons by protons . 

However, while the binding energies of nuclear 3- and 4- 
body systems have been theoretically computed and have led to 
some very important conclusions about nuclear forces (e.g. the 
forces operative between neutrons ), no very satisfactory 
scheme for mathematically treating collisions which involve 
more than two elementary particles has so far been given. In 
what follows some promising methods to solve such problems 
will be investigated and preliminary results presented. 

Theoretical calculations of collisions of light nuclei 
are worthy of interest for the following reasons: 

(1) They may confirm or weaken our confidence in existing 
nuclear theories; in particular, they may throw light on the 



1. H. A. Bo the and R. P. Bacher, Rev. M 03 . Phys. 8, 82 (1936). 

2. J. M. B. Kellogg, I. I. Rabi, N. P. Ramsey and J. R. 
Zacharia3, Phys. Rev. 57, 677 (1940). 

3. W. Rarita and J. Schwinger, Phys. Rev. 59, 436 (1941). 

4. W.. Rarita and R. D. Present, Phys. Rev. 51, 788 (1937). 

5. H. Margenau and D. T. Warren, Phys. Rev. 52, 790 (1937)* 

6. E. Gerjuoy and J. Schwinger, Phys. Rev. 61, 138 (1942). 

7. L. E. Hoisington, S. S. Share, G. Breit, Phy3. Rev. 56, 

884 (1939). 

8. D. Bohm and C. Richraan, Phys. Rev. 71, 567 (1947). 

9. R. D. Present, Phy3. Rev. 50, 635 (1936). 
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question whether nuclear forces are, partly, many partlole 
forces. 

(2) They may help to decide between certain ambiguities 
( such as the exact space and spin dependence of nuclear 
forces ), which the data at present available do not settle. 

Prom the very fact that more than two intimately inter- 
acting particles are involved it is clear that collisions of 
composite nuclei present a very complicated theoretical problem. 
While it is encouraging that the corresponding bound state 
problems have been solved with a reasonable amount of effort to 
within a few per cent, it may not be out of place to mention 
the additional complications of the scattering problems. The 
asymptotic form of the wave function, instead of having to 
vanish at infinity, must correctly describe the colliding 
particles. Above certain energy threshhold3 it must exhibit 
the possibility that part of the kinetic energy of collision 
may be used to raise one or both of the colliding particles to 
an excited state; or that an exchange of nucleon3 may take 
place, so that the resulting nuclei differ from the colliding 
ones; (e.g. V W-e* -t A 1 ) ; or, finally, that a triple 

or multiple disintegration may take place (e.g. La* 1 -4 P 1- —&■ 

( intermediate stages ) £ He.** * 1 )• Interactions 

with the radiation- and electron-neutrino fields will not be 
considered in this paper. We shall also exclude from our dis- 
cussion processes in which multiple disintegrations are 
energetically possible since in that case the asymptotic form 
of the wave function is very complicated and it was considered 
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advisable to concentrate first on a solution of the simpler 
problems • 



Experimental Data . 



In the last twenty-five years a great many experimental 
investigations of collisions of light nuclei have been carried 
out. Many of those performed before 1937 are discussed in 

lOjll 

Bethe’s and Livingston’s review articles on "Nuclear Dymanics n 
A few other typical references are given below . Scattering 

experiments with neutrons as of 1947 have been compiled by 

if 

H. H. Goldsmith et al. . Quite recently, collisions at very 

high energies have been investigated by means of the Berkeley 

lb 

cyclotron • 

It is our general aim to make possible a theoretical 
Interpretation of these experiments. However, this paper is 
devoted mainly to the development of theoretical methods and 
only very few numerical calculations have so far been carried 
out. Hence we shall not be concerned here with a detailed 
analysis of the experimental data. 



10. H. A. Bethe, Rev. Mos. Phys. 9, 69 (1937). 

11. M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 245 (1937). 

12. R. Schorr, J. M. Blair, H. R. Kratz, C. L. Blailoy and 

R. P. Taschek ( p-d scattering ), Phys. Rev. 72, 662 (1947). 

13. R. G. Nuckolls et al.( n-d scattering ), Phys. Rev. 70, 

805 (1946). 

14. N. P. Heydenburg and N. F. Ramsey, Phys. Rev. 60,42 (1941). 

15. H. H. Goldsmith, H. W. Ibsen, B. T. Feld, Rev. Mod. Phys. 19, 

259 (1947). 

16. J. Hadley et al.. Bull. A. P. S. 23, 15 (1948). 
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Previous Theoretical Work. 



Theoretical work on more complicated scattering problems 
has been mostly of a qualitative and phenomenological nature* 

It has been largely based on the resonance formula of Breit 

‘7 

and Wigner which exhibits the scattering process as the 
temporary formation of a compound nucleus with subsequent 
emission of the incident or some other particle. In terms of 
this theory it has been possible to understand the general 
features of nuclear scattering, but no detailed information 
about nuclear forces can be derived from it. 

Wheeler and Heisenberg have analyzed the general 
framework of collision processes without going into any detail 
of the interaction. Their theory serves two purposes: It 
elucidates the general features of scattering phenomena, 
regardless of the nature of the interaction, and it represents 
the first stage of any detailed calculation of composite 
collisions, based on some specific nuclear theory. 

Furthermore, as an approximation to the correct wave 
function, Wheeler has used functions of so-called group 
structure, which have been employed in all succeeding 
calculations. To determine the best possible functions of 
this sort, he made use of a variational principle. The 
contents of sections 1.6. and 1.7. are an extension of this 
method. 

17. G. Breit and E. Wigner, Phys. Rev. 49, 519 (1936). 

18. J. A. Wheeler, Phys. Rov. 52, 1107 (1937). 

19. W. Heisenberg, Zs. f. Phy3. 120, 513 (1942). 
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More directly related to the present paper are the works 

xo i\ to- i) 

of Schiff , Ochiai , Fltigge , Motz and Schwinger , 

2-t 

Buckingham and Massey , and Httcker who have calculated 
scattering cross sections in n-d and p-d collisions from 
approximate solutions of the- Schroedinger equation. 

These calculations are all based on a wave function 
with group structure, that is a function of the form 

(0.1) "Vi ^3 

where describe two collidors i and j when 

outside their range of interaction, and describes 

their relative motion. The sum is taken over all pairs of 
nuclei whose formation is compatible with the given energy 
and over all permutations of identical particles, with the 
signs adjusted in accordance with the exclusion principle. 

The functions T'*- are then chosen so that ^ satisfies the 
Schroedinger equation as well as possible. A function of this 
type has the correct symmetry properties and asymptotic 
behaviour; but in view of the strength of nuclear forces one 
cannot expect it to approximate too well the correct wave function 
in that region of the configuration space where the collidors 
interact, and it is just this region which determines the 
scattering cross section. 



20. L. I. Schiff, Phys. Rev. 52, 149 (1937). 

21. K. Ochiai, Phys. Rev. 52, 1221 (1937). 

22. S. Fltigge, Zs. f. Phys. 108, 545 (1938). 

23. L. Motz and J. Schwinger, Phys. Rev. 58, 26 (1940). 

24. R. A. Buckingham and H. S. W. Massey, P. R. S., 123, 

179 (1941). 

25. R. A. Buckingham and H. S. W. Massey, Phys. Rev. 73,260 (1948). 

26. K. H. HOoker, Phys, Zs., 43, 226 (1942). 
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Of the calculation mentioned above, those of Buckingham 
1H if 

and Massey ’ are the most extensive. Their method is in 
essence the following: 

The centre of gravity motion of the two neutrons and 
proton, participating in an n-d collision, is separated out 
and there remains the following equation for the wave function 
in the relative coordinates: 

(0.2) ^ Sk' , 



Here T is the kinetic energy operator for the relative motion 
of the three particles; i3 the interaction potential 

between particles i and supposed independent of the position 
of the third particle; is the energy of relative motion of 

the colliding neutron with respect to the deuteron in the system 
in which the centre of gravity of all three particles is at rest; 
and is the binding energy of the deuteron. 

Let us introduce relative coordinates to which we shall 
have frequent reference in the following pages. If 
and are the absolute position vectors of the two neutrons 
and the proton, then P v , ^x. » are defined by 



( 0 . 3 ) 






and similarly f. 



^ are defined by 
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(0,4) 

’ t - 



r, - 



su 



2 . 

— ^jjJL* 

X. 

R x + fti. 




In terms of these coordinates T has the form 

(0 ‘ s) T-- LV + i %) 

or the two other forms obtained by replacing 1 by 2 or 3* 

M is the mass of a nucleon* 

The wave function is a function of the six relative 

coordinates represented by ^ and ^ ( or any two linear 
combinations of these ) and of the spins of the three particles 
Equation (0.2) is to be solved subject to the condition 
that when either or become large, ¥ describes a 

deuteron and a free neutron colliding with it. Thus 

^ f (?») (e 1 * 1 '' 1 + ^ ' ,;W ' 






Oji f -$> 



cu r. 



oO 



£?0 . 
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where is the ground-state wave function of the deuteron, 

t^lkl»C ; , f is the 

scattering amplitude and and ^^are appropriate spin 
functions which differ according to whether we consider quartet 
or doublet scattering. 

Although we are not familiar with a rigourous proof that 
a six-dimensional partial differential equation of elliptic 
type has a solution which is finite everywhere and satisfies 
such asymptotic boundary conditions, one intuitively expects 
this to be the case. However a reasonably rigourous solution 
is undoubtedly very complicated. Buckingham and Massey, to 
simplify their problem therefore assume that the exact solution 
can be everywhere well approximated by a function of the form 

(0,7) £ = 

This i3 equivalent to allowing only for the distortion of the 
free neutron function ( since P is as yet undetermined ), but 
disregarding the distortion of the deuteron by the neutron. 

To find which function F makes ^ the best possible 
solution, they substitute (0.7) into (0.2) and, using the wave 
equation for the deuteron, they find that P should satisfy, in 
the case of doublet scattering, 

* V)\ V (r.) 

-- (^ - -vV« £f(f # ) 



( 0 . 8 ) 
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Here w, h, m and b represent the strengths of the Wigner, 
Heisenberg, Majorana and Bartlett forces# 

For any fixed value of ^ , (0*8) could in principle be 

solved for FM giving' ^(r). But as there is no reason why F 
should be independent of ^ it is impossible to satisfy (0.8) 
by a single function f'(r) . This is just an expression of the 
fact that a function of the type (0.7) cannot solve the wave 
equation exactly. 

To obtain, as it were, an average solution for F, 
Buckingham and Massey multiply (0.8) by and integrate 

over the ^ -space. This leads to an integro-differential 
equation for F, of the form 

(0.9) ^ 

+ (au- -A) j Q(r,r') F (r-)[Ur'l 

where U and Q involve the known interaction potentials and 
the deuteron function. This equation is finally solved by 
numerical methods. 

Previous to Buckingham and Massey, Motz and Schwinger 
have calculated the scattering of neutrons by deuterons in 
the zero energy limit. The essential features of their 
calculation are similar to those of the British workers. 

Thus, they too consider the deuteron as rigid and are led 
to an integro-differential equation for the function F. 

To solve it they replace it by an integral equation according 
to the scheme: 
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(y\a v J> F(r) s r "3 ^( r " ^ 



(0.10) 

F (0 



A k- r 

e -t- 



HTT l r-Vl 



After carrying out angular integrations analytically the 

remaining integral equation is solved approximately by 

replacing it by a set of simultaneous linear equations. 

At T” | \ 

Finally, Iiticker has tried to write r|rj as a power 
series inland determine the coefficients to best fit the 
Schroedinger equation. However his results depend, by a 
factor of about 1. 5-2.0 on the point where the power-series 
is taken to pass over into the asymptotic form of F. 

The other calculations do not contain any essentially 
different methods. 

The different authors using various kinds of ( central ) 
interactions obtained widely different results. For the 
total cross section of n-d scattering at 0-energy they are in 

XO 2.1 

units of barns as follows: Schiff : 6.11 ; Ochiai : 3.2 ; 

Motz and Schwinger : 4.57, 6.91 ; Buckingham and Massey : 

16 

3.0, 3.9 ( extrapolated ) ; HScker : 1.88, 2.95 . The 
experimental value is near 3.3 ) 

Other two- and three-particle systems do not exhibit such 
a strong dependence on the type of interaction, at low energies. 
If taken at their face value, then, these results would • 
provide a useful method for distinguishing between various 
kinds of interaction, all of which are compatible with other 



27. L. J. Rainwater et al., Phys. Rev. 73, 733 (1948) 
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low energy phenomena. 

However, such a procedure would be without firm foundation 
in vlev/ of the unreliability of the group structure function, 
on which all calculations are based. Thi3 emphasises the 
need for a systematic method of solution, which is not 
restricted to any particular kind of function. 



Summary of the Present Paper . 

. It is not the primary purpose of this paper to add 

another approximate calculation to those already published, 

but rather to investigate the possibility of a systematic 

procedure for solving the Schroedinger equations corresponding 

to composite scattering processes. It was natural to look 

for a variational method, as this has been the chief tool 

for solving the problems of the bound 3- and 4-particle 

systems 1 and has recently been applied by Schwinger 

to 2-particle scattering problems. 

There are two types of practical variation procedures: 

3 * 

One is the Rayleigh-Ritz method of undetermined coefficients; 
the other the variation iteration method employing Green's functions. 
We have tried both methods in the many particle scattering 
problem, and we strongly incline to the opinion that the 

28. N. Svartholm, Thesis, Lund, Hakan Ohlssons, Boktryckeri (1945). 

29. Unpublished lectures. 

30. w. Ritz, J. relne angew. Math. 135 (1909). 
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Rayleigh-Ritz method la the more useful. 

Using this method we have made exploratory calculations 
on 2- and 3-particle scattering problems, some with the 
inclusion of tensor forces. On the basis of these we believe 
that this procedure can indeed be used for a systematic 
investigation of collisions of light nuclei. 
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I# THE RAYLEIGK-RITZ METHOD IN COLLISION PROBLEMS. 



1.1. General Remarks. 



In his original paper on wave mechanics, Schroedinger 
derived the wave equation from a variational principle 
analogous to Hamilton’s Principle in classical mechanics. 

The variational formulation of the equations of motion of 
wave-fields has been very fruitful not only in theoretical 

n 

developments of field theory but also in approximate 

33 

solutions of complicated quantum mechanical problems . 

The method most extensively used in practical problems 

jo 

has been the Rayleigh-Ri tz method of linear trial functions. 
However its application has been limited to systems in bound 
states, whose wave functions vanish rapidly at infinity. 

Recently, Heisenberg has pointed out the correspondence 
between the phase shift energy relation of a collision problem 
and the energy spectrum obtained if the participating particles 
are enclosed in a large box. 

Consider, as the simplest example, a particle of mass M 
under the action of a 3hort range scattering centre. The wave 
equation for the radial part of the spherically symmetric 

I ' 

solution Is 

31. E. Schroedinger, Ann. d. Phys. 79, 361 (1926). 

32. See, for example,. G. V/entzel, "Einfilhrung in die 
Quantentheorie 3 . Wellenfelder", Wien, Franz Deuticke 
(1943). 

33. E. Hylleraas, Zs. f. Phys. 54, 347 (1929). 
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(i.D ( = €uW ; 

l x M ^ J 

where U (a^ la the scattering potential and E the total energy. 
Asymptotically the solution becomes 

d.2) tc — ^ A + 

where the phase shift \[ e ) , corresponding to the energy E 
of the incident particle must be found by a detailed 
solution of (1.1). 

Now let us ask for the energy levels of the same 
particle in the same potential field when enclosed in a 
spherical box of large radius, a, on which the wave function 
must vanish. Since the differential equation is unchanged, 
all those solutions of the scattering problem will satisfy 
the bound state problem for which 




It is clear that if we have solved the scattering 
problem and determined the function 'll 6 ) we can find the 
energy levels, of* the bound particle by solving (1.3) 

for Ej and conversely if the spectrum, , of the bound 
particle has been determined the phase shifts ^ (€x) can 
be found from (1.3). 

These considerations lead us to expect that the 
Rayleigh-Ritz method, which has been so valuable for 
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calculations of energy levels of bound systems, will apply, 
with certain modifications to the calculation of phase shifts 
and scattering amplitudes. 



1.2. Three Related Treatments of One Dimensional Problems . 
The Collision Determinant. Numerical Illustrations. 



In thi3 section v/e shall illustrate our methods by 
treating the scattering of a particle by a fixed, short 
range scattering centre. We shall show that thi 3 problem 
can be attacked In three closely related w ays. 

Let us first simplify the Schroedinger equation (1.1) 
by introducing a convenient unit of length, , and writing 





x. 

* 



It then becomes 

( 2 - 2 ) u (*)-0 



It will be our problem to solve (1.5) subject to the conditions 
that 

(2.3) 0) « O j U.(x) 3 StUx VC* -Tfou. ^ ^ * * 

Od — ^ o'? . 

In particular we are interested In the phase shift which 
determines the scattering cross section. 
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Method 1: The Energy as Proper Value . 

We know that (2.2) and (2.5) have a solution for every 
positive « However, instead of looking for the phase shift 
corresponding to a given K , we can reverse the problem and 
ask for the K corresponding to a given logarithmic 
derivative, L, of u at some large distance a. This means we 
look for the proper value belonging to the differential 

equation (2.2) and the boundary conditions 

(2.4) U.(o) ^0 j l\' [ A) ^(<0 . 

This is a proper value problem of a standard type and it 

is easily verified that as calculated from the equation 

( 2 * 5 ) ^ ^ - U* X -fc- 14. w V U?" | L U - [<xj O o 

is stationary with respect to all variations of u which 
satisfy the condition <5u(ojcO. 

The Rayleigh-Ritz method is directly applicable to (2.5). 
We simply try for u a linear combination of trial functions 

(2.6) W, Hr - - — 

( with all U.^ (/?) -0 ) and substitute this expression into 

(2.5) * This leads to an equation of the form . 



(2.7) Q C ** C ) ~ 0 



) “ » 



where A;\ . 6 are integrals involving CtwcL tc\ 

U 2j 



34. R* Couront and D. Hilbert, Methodon d. Math. Physik, 
vol. II, 508, Berlin, Spi'ingor, 1937. 
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The stationary property of K with respect to variations 
of the leads to the following system of equations for the 
determination of these coefficients: 




Ast l 



jij - - • n. 



As in the case of the bound state problem this system is 

It 2, 

compatible only if K has certain proper values which 

can be found by the solution of the corresponding secular 
equation. They are approximations to the correct proper values 
of (2.2) under the boundary conditions (2.4). Their errors 
are of the order of the square of the errors of the trial 
functions . 

Now the phase shifts of the exact solutions corresponding 
to the given L satisfy the equation 

(2.9) k -t -L . 

If we substitute our approximations in this equation we 

obtain, to the same ( namely the second ) order of accuracy, 
values for ^ ^ corresponding to the energies \ijj~ 

Then, either by varying the pre-assigned value of L or the 
,, cut-off ,, distance a, we can obtain for all other energies. 

In this form of the variational method, we have followed 
very closely the analogy between bound and unbound systems 
pointed out by Heisenberg. However the quantization of the 
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energy which entails the necessity for solving a secular 
equation of the usual kind seems artificial and can indeed 
be avoided* 

Method 2« The Logarithmic Dorlvative a3 Extremum . 

Since for a given L, K. 1 ” as calculated from (2.5) is 

stationary, so also mu3t L be stationary when calculated with 
z_ 

a given K . Using the trial functions (2.6) we are led 
again to (2.7), but now we take as fixed and L as 

adjustable so as to make the ensuing system of equations 
compatible* Therefore we exhibit L explicitely and write 
Q in the form w 

(2.10) 

~ — 0 

It is convenient to transform the coefficients c • to a new 

A 

set d * as follows: 

A* 

^ l - bv. v C 6 ') ^ s C 

( 2 . 11 ) 

^ ~ c ^ *1 - 5-) 3, . - n 

In terms of the d ; , Q has the form 

A 

(2,12) Q - 2 D- ^ j -t u <£ ^ - o 

-A - \ 




The stationary property of L with respect to variations of the 
d \ now leads, to the equations 
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^ [vO') -i L_ d v ^ 

(2.13) * 

P;j U v ) d'3 — o J ^ *. 

j' 1 

To make the system of equations (2.13) compatible L must 
satisfy the dete.rminantal equation 



(2.14) A £ 



L*v 0 U 




I 

I 




Du- - 

Div. 



1 

I 

4 

1 





\ 

i 

\ 

I 

) 




0 . 



• • 

Hence L is a quotient of two determinants, each of which is 

a function of K • 

Let us note that A is symmetrical and of degree 1 
in L, in contra-distinction to the secular determinants in prob- 
lems of the bound state, whose, degree in the unknown energy 
is as high as the number of trial functions. We shall see. 
later ( sections 1.5. - 1.7. ) that the structure of A. or 
of the corresponding determinant occurring in method 3 reflects 
the type of the collision and for thi 3 reason call it a 

3 r 

colli sion determinant • Its dogree in L Is of course 

35* It is related to Wheeler 1 a Fredholm do terminant , ref # 

18, eq« (45) ♦ 
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Invariant under the transformation (2.11) which was introduced 
only for reasons of greater clarity« Indeed one can easily 
verify that the determinantal equation for L, as deduced 
directly from (2.10), 

(2.15) | 4 >- “•> l a ) \ -° 

is also of the first degree in L, since the matrix 
\\ u 'j (<0 ll of the coefficients of L has rank 1* 

There remains as an unsatisfactory feature of this 
treatment the dependence of the phase shift, deduced from 
and L, on the cut-off distance a, even when a is 
well beyond the range of interaction. This question will 
be discussed in more detail further on in this section. 
Meanwhile we shall indicate another method which doe3 not 
have this drawback. 



Method 5. The Tangent of the Phase Shift a3 Extremum . 

Let U3 return to equation (2.2) and ( writing £ 
instead of ) note that, in the bound state, £ 

calculated from 



(2ol6) X 




o 

is stationary relative to all variations of u with 
and Su(©o) :=• 0 • For, 



<5i*(o) =o 




(2.17) 
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and both the integrand and integrated term3 vanish. 

In a collision, however, the wave function has the form 
(2.3) at infinity# If we restrict our trial functions to 
those having the correct asymptotic form, although . ' . 
possibly v/rong phase shift, then 



(2©18) S U — & 'AaX~\ v| , ev-o K x 



X — ^ oO 



Hence 



d (I 



(2#19) 


m. "T u. 

oU 






7L* oO 


so that 


S[X + 




(2.20) 


K yfccw u J - 0. 



14 d 




Since for the correct wave function I r. 0 , it follows 

that the equation 



(2#21) 3T *V" 1C vj^ — ^ Xklt*. Wy 

where is the trial phase shift, is a stationary 

expression for tan v| • 

To apply the Rayleigh- Ritz method it is essential to 
write this variational principle in a homogeneous form# 

Let us therefore consider wave functions of the asymptotic 
form 



36© L# Hulthen ( Extrait, Dixieme Congres des Mathematicians 
Scandinavos, Copenhague, 1946 ) has used an equivalent 
equation to find >| by further restricting his v/avo 
functions to satisfy the condition I z. 0© 
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(2.22) U ? ft Rx -t 6 i^>V.oc 

instead of (2.3). We then find in the same way as above 

(2.23) 14 A B )<- 

as a homogeneous variational principle for tan V| . If we use 
a trial function of the form 

(2.24) H, * A Sw 4 C L k l+ C K k 

where 

(2.25) k L 4 c 3 -t ' “ — B >4)C 

• ^ 

then since tan ^ multiplies only A , substitution of (2.24) 
into (2.23) will give a quadratic form of' the same structure 
a3 (2.12) which in turn will again lead to a linear collision 
determinant for tan i| . 

The drawback of this formulation is that the trial 
functions must be of the form (2.24), (2.25). For practical 
purposes it may therefore be preferable to re-adopt a finite 
cut-off as follows: Let a be a distance beyond which the 

interaction i3 negligible and let us use as trial function 

C | Uk l" fC L U 2.-*' C \^v\ s X < Q. 

(2.26) 

^ _ l\ C ° * - ft KX. + ft VCX 

where 

ft ■= | ( a) u a + d l) ( q) • ua j 

6 , ( ~ V uh>l ( a ) 5,v 



(2.27) 
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A and B are defined so that u is smooth across x Clearly, 

1 °) 

since u has the form w- for x > a, the integrand of I is zero 
for x a ( see (2.16) ) so that 



(2.28) 





When (2.27) and (2.28) are substituted in (2.23), we have a 
variational principle involving only [x) , whose asymptotic 

form no longer concerns us. The price paid, of course, is the 
reappearance of the cut-off distance a. 



The Well-Depth as Proper Value . » 

Before going on to a discussion of these methods let us 
note that if we write V/ (x) ^ -^(x) • where f(x) determines 

the shape and ^ the depth of the potential well, we can 
regard ^ a3 the proper value in the Schroedinger equation 

(2.29) J + VC” -V ^ k W = ° • 

l fa' J 

^ ; as calculated from the equation, 

°° /V 

(2.30) \ \< V ^ ^ ^ 0 

is stationary if 

# 

(2.31) du(oj r. 0 ~ O . 



This idea was used in connection with bound state problems 
by Svartholm ( ref. 28 ) . 
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Both K and are now regarded as known while the strength of 
the interaction is the unknown. This is in contrast to 
methods 1-3 in which the interaction was regarded as given 
and U and or K and L were determined as functions of one 
another. 

While this approach has no advantages in connection 
with the Rayleigh-Ritz method, it provides a useful view- 
point in the variation iteration method to be discussed 
later ( section II.8.). 

Connection between Methods 1-3 : 

Mothod3 1 and 2 differ only w ith respect to which of 



the quantities L and K is presumed known and which is 
unknown. With the same trial function and the same a they 



On the other hand method 3 differs slightly from 1 and 
2. Let us trace the connection between the cut-off independent 
form of method 3 and method 2: 

In the latter, L is determined from (2.5), which after 
an integration by parts gives 



where I is defined by(2.28). If u has the asymptotic form 
(2.22)wo can evaluate the two boundary terms, which leads to 



must lead to the same relation between tan and K. . 



(2.32) - U ^ 
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(2.33) 



A „ / 

S‘/vt^K^ + wj3(J ( K. <\ 4 *1 * ) -+IE 

C^> V|^ l \ ' 

o'- 

+ R CVt" (_ VC 0. -4 V|\ 5 W 1, fx*-* n*) - 0 

Ov’V ' 1 J *. 

I c~-> y\r 



) 



where ^-'ta--' ( &/*) • Kultiplying (2.33) by 
we obtain 



(2.34) 



SU^Ka-^*) 

- o 

Su l [ Kd + hi) 

At the particular cut-off distances where 1 l 

(2.35) Ka - 

a / 



(2.34) becomes 



(2.36) 31 4 A X !/t<Xv> v^, — A K txbv^. 



fhis agrees precisely with the result (2.23) of method 3. 

Y/e now ask how ^ as determined by method 2 varies with 
a for a given trial function u which ha3 the correct asymptotic 
form but incorrect phase shift. In (2.34), let us use the 
abbreviations 



(2.37) 



^ s K* + } Svj r I'j * - ^ 

\ "X , t. 



By some elementary manipulations we can then reduce it to 

(2.38) UX <5 h - cvl^-'/O, 



For the correct wave function, J -s- 0, so that by (2.38) 



I. Co$ ^ — 
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d O ^ except where cot z =. + ©o , that i3 where 

^ * At those points S is undefined. This 

is understandable from (2.5) since at those points u = 0 so 
that L is not determined. 

The variation of with z is shown for different 

values of J in Graph 1. For 3mall values of J, cot z is 
negligible as compared to l/J, except at the danger points 
UA-4^ ~ ^ ~ * 30 ^at in ^he ii m it as J — ->0, & ^ 

becomes independent of a. 

Bilinear Formulation. 



The Schroedinger variational principle for bound systems 
is usually stated in either of the two forms: 

i I'm) - v (*) ■+ <**■-<> 



(2.39) 



JCUO 



or 



(2.40) 



cO 

Q = \ u l i - \j[y) -v <c \ u ^ 
O l d** J 



where all trial functions satisfy 
(2.41) U [o)=.0 j 5u[ot>^ = 0 - 



Let us recall that a modification of (2.39) for the 
case of a collision led to methods 1 and 2, while (2.40) 
was the starting point for method 3. 
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We can replace (2.39) and (2.40) by the following 
equivalent statements, involving bilinear instead of 
quadratic forms: 





a* * H-I 


[ t )l 


f i 


^ _ v/(*) -t ^ ^ fix ^0 


(2.42) 


J o l 


\&07l ) ( 


s. d-^c 


or 


<5 Q a ^ o 

6 ix ~ ^ k , ^ 




~ \/(x) + k'b ^ » 0 


(2.43) 








where 


5 Q 12 . - o 

all trial functions U.^ 


and 


satisfy 


(2.44) 


=: 0 ; 




lr.0 j / s ^2. . 


S in 


(2.42) and (2.43) 


indicates 


independent variations of 



and U^. The truth of (2.42) and (2.43) is directly 
verified from the Schroedinger equation and the boundary 
conditions (2.44). 

Clearly the corresponding - variational principles for 
collision problems can be similarly modified. 'Without 
going into details we shall merely state the results. 



Methods 1 and 2: 



(2.45) 



i\- (.*)(%) 

-t L* u, (a) ux ( $) - o 



u, 








3 vipj, ~ 0 

if the trial functions satisfy 
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(2.46) ^;| 3 ) ;0 ^ X o ~L . 

( Compare to (2.5).) 



Method 3: 

,2.47) H lt S It, =° 

Iii = ^ u I | H R 1 ' - ^ (.*) j Uj, <=£*. 

(2.48) £1^-0 



if the trial functions satisfy 

3 wu(o)=o ( 

( 2.49 ) ' 

tC ^ /^x» S i.'-*'- vt X- — t 6*^' c-->-5 V4 of 

( Compare to (2.23), (2.24).) It is not at once apparent 
that tan as calculated from (2.47) is symmetrical In t4. ^ 
and l* • That this is so, however, follows from the fact that 

oO 

^li = ^ u ijj'v + K v " V(*)| u i 



( 2 . 50 ) -- [ 4 rf- Vl»)\ •*, ^ + 






id 



In one dimensional problems there is no advantage in 
using the Raylelgh-Ritz method in conjunction with the 
bilinear variational principles instead of the quadratic 
one3. But wo shall find the bilinear version useful in three 
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dimensional problems ( sections 1.4. - 1.5. ), where the 
Schroedinger equation has more than one solution, and in 
connection with the variation iteration method ( sections 

II. 8. - II. 10. ) . 

Numerical Illustrations . 

Using the methods 2 and 3 we have calculated 3 S 
scattering of neutrons by protons at zero energy. We have 
assumed the interaction 

\/(x) := -^.01435* 1 

(2.51) 

- 0 X >| 

( x measured in units of /i 0 = Z-30 jtio”'* cm )> v/hich give3 
the correct binding energy for the deuteron. The wave 
equation, in this case, is simply 

it 

(2.52) 23 _j. y/( a ) k. = o. 

cU: 1 ' 

and can, of course, be solved exactly. For x > 1 the 
solution has the form u - A ( x + X ) where X is the desired 
limiting ratio of tan vj / £ , as K approaches zero. 

As trial functions we have used 

h 

(2.53) -s *2. 

I 

with n z. 1, 2 and 3. Since (2.53) does not have the correct 
asymptotic form. It was necessary to use method 3 with a cut- 
off, a. The interaction vanishes at x * 1, which is the most 
appropriate choice for a. However to study the dependence of 
X on the cut-off distance, vie have made calculations also for 
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larger values of a. Clearly a function of the form (2.53) 
can approximate the correct wave function better in a small 
than in a large interval, so that as a increases we expect 
our approximation for X to get worse. 

The results are listed in Table Is 

Table I . 

Approximations to X s / K < 



X a - 2.080 



n 




c i 


X 




CjX H 


K c^x*" 


C X + c 

1 1- 


x- , 3 

X + C X 
3 


u 


Me thod 2 


Method 3 


Method 2 


Method 3 


Method 2 


Method 3 


i 


- 


3.957 


+ 


1.338 


- 2.112 


- 2.126 


- 2.083 


- 2.084 


2 


+ 


4.043 




1.338 


- 2.858 


- 4.039 


- 2.201 


- 2.206 


3 


+ 


2.416 


+ 


1.338 


- 6.593 


- 51.606 


- 2.385 


- 2.520 


4 


+ 


2.011 


+ 


1.338 


+968.870 


+ 7.634 


- 2.936 


- 3.315 



These results exhibit a number of significant features: 

1. The methods 2 and 3 lead to errors of the same 
order of magnitude. ( In the present example method 2 is 
somewhat better. ) 

2. A3 expected, the results grow progressively worse 
with increasing a. 

3. The dependence on a becomes weaker the larger the 
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number of parameters ( except for the 1.338 column; see 5. 
below ) . 

4. There appear singularities, owing to the wrong 
asymptotic behaviour of the trial functions. As the number 
of parameters increases the singularities move outwards, to 
infinity. 

5. The trial function c ( x, in method 3, leads to a 
result independent of a but of the wrong sign. This can be 
understood as follows: c^ x is the correct solution corresponding 
to no interaction. It has ( accidentally ) the correct 
asymptotic form ( though wrong X ) and therefore, by the theory 
of method 3, must lead to results independent of a. In fact 

we have ^ 

T - j - V(x) oloc 

(2.54) 0 , 

ft c 0 

all. of which are independent of a. Using these in (2.23) we 

obtain for the determination of X - tan y] / K, 

Jo 1 

(2.55) - ^ (*) etc C.^ 

which is precisely the Born approximation for tan ^ . For 
small V the latter gives good results. But as V increases 
( tan V| has a fixed sign while the true tan 

eventually changes sign. This has occurred in our case, 
hence the v/rong 3 ign in the 1.338 - column. 

In summary we may say: 

(a) Method 1 is the least natural to collision problems. 
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(b) Method 2 appears to give slightly better results 
than method 3; but it suffers the aesthetic flaw of necessarily 
involving a cut-off distance. 

(c) Method 3 is closely related to method 2. It does 
not necessarily involve a cut-off and yields the Born 
approximation for the phase shift as a special case. 

(d) The Rayleigh-Ritz principle applied to either 
method 2 or 3 leads to a collision determinant instead of the 
usual secular determinant. 

1.3. Variation of the Phase Shift with Energy. A Stationary 

Expression for the Width of a Resonance Level . 

7/e have seen that methods 2 and 3 lead to stationary 
expressions for the phase shifts corresponding to arbitrary 
energies* It follows that it must also be possible to 



the variation of the scattering cross section with energy. 

It Is indifferent which of the tv/o methods we use, but 



obtain a stationary expression for 




which determines 



method 2 seems more convenient. We introduce the definition 



(3.1) 





-ti 



where is the logarithmic derivative of the correct 

function belonging to the energy K**. By the Schroedinger 



equation 




and by the stationary property. 
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(2.5), we have rigouroualy 

(3.2) ~ s oOh), 

for any function tc •=. u. ^ S u. ^ . Differentiating (3.2) with 
respect to K* we obtain, to the first order, in S^k , 

cL ( L\} 

(3.3) ___JLL -L - o 

0 (. l< 

Hence by (3.1), we find the equation 

(3.4) 2K \ ^ 4 — p 1 K ( a ) “ 0 

) d K 

o 

which is a stationary expression for j ck K 

From this equation we can determine as follows. 

We note that 

(3.5) L |4 ■=• & (K.d4^w.) 



Hence 

, , c ’ Lk ( V. o>+ M v.) - -+K. 

(3>6) ITU - v 1 1 A d«J 



so that 



(3.7) 



a* “Ttil 

iAL-k. S\^f uq-f + __ ei 



sLv. 

<K 



4 K.’ 






^ U V dx 



- « 



14 u v («) 



, S*/k” (/a -4 



4k u 

a 

*U< 
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We may observe that this stationary expression for 
involves only the wave function and the correct phase shift 
, while all reference to the interaction has disappeared. 
Let U3 check that if ^jc/^fcjis computed from (3.7), 
using the correct wave function , the apparent dependence 
of ^‘|K.[^ v J°n a disappears, as it must, as long as a exceeds 
the range of interaction, . 

For x > & 0 , the correct has the form 

( 3 . 8 ) C 



Hence 

cL i| n. 

< 3 - 9) 7i ?) 




H- tt 1 - 




<2 K. 



Differentiating with respect to a we find 



1.1 

V d(*y)J 

(3.10) 



\ 



KC* 



. C 





-V^k)) Hr ^ S ( R a -t 



0 

which proves that d(n v ) I s independent of a. 

As an application of (3.7) we shall derive a stationary 
expression for the level width at resonanoe. 

In the neighbourhood of a resonance energy, defined by 
the condition that x > tan has a singularity 



of the form 
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(3.11) 





where i 3 the resonance energy. C 1 13 called the level 

width because the cross section in the vicinity of i3 

given by 



(3.12) cr ^ 



Nov/ near K.^ ? 
(3.13) 
Therefore 





(3.14) ~Uw ~ 




Hence, comparing (3.11) and (3.14) and using (3.9) with 
we find the following stationary expression for ^ T" 1 j 





(3.15) 






^ *]k. 

*(•£.) 






. 4k 



t ^ / 

^ J tv oUc -f K.) 
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1.4. A Variational Principle for the Scattering Amplitude . 

'When dealing with a wave function of given angular 
momentum, a correspondence to a bound state problem can be 
established because we know, a priori, that asymptotically 
the solution has a radial logarithmic derivative independent 
of the angle. But If we want to treat the whole incident 
plane wave at once no such information Is available. We 
have already seen that only method 3 need not involve an 
enclosure and that the Born approximation for individual 
phases can be deduced from it. Both facts suggest that It 
v/ould be the most promising attack on the scattering amplitudes 
of three dimensional problems. 

The wave equation for the relative motion of two colliding 
particles of equal mas 3 M Is 



(4.1) U(k)j Sf Ik) - e 



where R is their relative position vector and U (R) is the 
interaction potential. Setting 



(4.2) 





where is our unit of length, it becomes 
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(4.3) ^ -* K* - V[*)\ 't(t) =»0 

We shall nov/ follow a procedure corresponding to the bilinear 



form of method 3 ( see (2.47) - (2.49) ). We define , and 

as the two solutions of (4.3) which correspond to plane 
waves incident along -K l and \< u and have the asymptotic forms 

„ A l< A 

(4.4) 

Here [ - [{^1 and is the scattering amplitude 

in the direction K of a wave incident along . Corresponding 

to (2.47) we now define 



4 / -k 1 - 


Xvt.. r 

$> e -v 




* *Vr 

— $> £ + 



A 

X 



-o 



(4.5) 1, 

K 1 



[y -t k\- V^) (Jr) 



In virtue of (4.3), we find after an Integration by parts, 

<4 ‘ 6) SI K h «, \ +-K,1 

s 

£ 

where S is a large sphere, and -- denotes differentiation 
along the outward normal. 

Now let us admit only trial functions of the asymptotic 
form (4.4), but with possibly Incorrect f T s, so that 
asymptotically 



(4.7) 5^(k>,, *) 



KA 



-A 
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where <5 (. K.) is the error in the scattering 

amplitude of the trial function Hr**- Clearly only the plane 
wave part of con trihutes to the surface integral (4.6) 

so that 



SI, 









a vL A 



(4.8) 



-t K A ^ 

* -i 4L 

s\ 2 w 



iS. 



M C o a >, \ 

Dirac has shown that the factor multiplying oA ft^Kjhas 

a & - function like character for large values of /i , 30 

that the integration gives simply 



(4.9) -- 4 l» S (Ki , l<i). 



Since for the correct 4X’ -L — 0 , it follows 

that (K, *0 > as calculated from the equation 



(4.10) 



X 



A*) 

^2^1 + 4 tr \ 



C -=^ 4"rr *0 



is stationary with respect to independent variations of 
and , subject to the conditions (4.7). 

Symmetry Properties . 

Equation (4.10) involves two trial function, and 

• These differ by a rotation in the coordinate space 

37. P. A. M. Dirac, The Principles of Quantum Mechanics, 

( Clarendon Press, Oxford 1947 ), third ed., p. 191. 
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and, in general, also in their functional form ( to the 

first order ). Let us indicate the functional form Ly super- 

, 10 10 

scripts and write explicitly ^ for our trial 

functions, and )] (lA L ^i4 v Jfor the f obtained from (4.10) by 
(A 3 C* } 

using first and second in (4.5). Integrations 

by part3, analogous to (2.50), reveal the following symmetry 
properties which are not at once apparent from (4.10): 

(4.11) K. t ) - , 

(4.18) (-**,*») ^ { ''"C-**')' 4 -)- 

(4.11) implies that f is Invariant under interchange of the 
functional forms of , and * as be verified 

below in the case of the Born approximations. (4.12) 
guarantees that the scattering kernel \<^ Is symmetrical 

when computed by (4.10). 



Connection with Born Approximations . 

The simplest admissible trial functions are 



(4.13) 



't-Ki * e. 



-X kj . r 



Mr ^ e. 



x K t . r 



38* Ref. 18| oq. (2.52) 
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Substituting these in (4.10) we obtain 

(4.14) ^ ~ Mw 

which 1 s identical with the result of the first Born 
appr oxima ti on • 

We obtain the second Born approximation if we use one 
first Born approximation trial function and one plane wave; 



i.e 


. either 










IM 


-a 1C. V“ 


(4. 


15) 1 1 


* Hr-*. 5 


e 




Mr 


o; 


a yu*v 
e ^ . 


or 












12.) 






. mt- 


H, - +-*, 




(4. 


16) 










(t) 






H' 






in 


agreement 


with (4.11). 


Higher 


obtained in 


a similar way. 








Thus, apart from re-deducing the Born approximations 
from a variational principle, nothing new is gained by using these 
iterated functions in (4.10). 



Application of the Raylelgh-Rltz Method . 

We can use the Rayleigh-Ri tz Method in conjunction 
with (4.10) to determine the scattering amplitude 
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This is analogous to method 3 ( section 1.2.) for calculating 
the tangent of the phase shift. For simplicity, we use the 
same functional form for ^ an< i and write 



(4.17) 



k,. ** 

t-K, * C., e +c t u. Vl -e 



i 1 '/ 



^ '<>, - c i «- 



K.,. r 



+ c i j a, -f . . .. U. , w 






where the trial, functions U have the asymptotic form 



(4.18) U.(k') 






^ {*> «') 



: v. a. 



Clearly for use with trial functions of the form (4.17) instead 
of (4.4), (4.10) must he replaced by the equivalent homogeneous 
equation 



(4.19) 1^,4, 1 C; {4*,,*,} «45TC ^ 

where is defined by (4.5). 

When (4.19) is written out in terms of the trial functions 
it assumes the form 

U 



(4.20) “E C -i C ) =- M« K,) C, 









where we may take j =- • The stationary property of f 

relative to variations of the C* then leads to the equations 



(4.21) 



f m , A 'j c i - ^ {l *,)■*, 

K 

ft;j c j - P , t K 



1-1 
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They are compatible only if the determinant of the coefficients 
vanishes, that is if 

Aji, - - - Am 

(\ lx A it1 



(4.22) = 



A K » 



I 

Akv » - - - - Auu 



- 0 



f( , VCj ) is uniquely defined by (4.22). The fact that it is, 
in general, complex is reflected by the collision determinant 
w ^ ose matrix is symmetrical but not hermitian. 

On varying the angle between and one can obtain 
stationary expressions for the scattering amplitudes, and 
hence for the differential cross section, in all directions. 



1.5. Mon-Central Forces . 

It is well known that if two particles interact by means 
of a potential which depends only on their distance of 
separation, the relative angular momentum is conserved. This 
fact is the basis of the phase method in collision problems, 
y/hose variational formulation has been discussed in sections 
1.2. and 1.3. 

However, the hypothesis that nuclear forces are purely 
central is untenable in view of the quadrupole moment of the 
deuteron, and has been replaced by the assumption that they 
are a mixture of central and non-central forces. Rarita and 
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Schwinger have worked out the combinatorial aspects 
introduced into the 2-body problem by the existence of such 
non-central, so-called tensor forces. In particular they 
have shown that the angular momentum is no longer conserved, 
and that the eigenstates of the Hamiltonian, instead of being 
S, P, D, . . . . states are now, in general, mixtures of such 
states; Most important, in the neutron-proton system at low 
energies, is the S + D mixture. 

In our method of section 1.4. , the conservation of 
angular momentum did not have to be assumed and thus it 
applies v/ithout change to collisions involving tensor forces. 

This method deal3 with all angular momenta at once. It may 
hov/ever be desirable to treat each eigenstate separately. 

This can be done in a way related both to the three dimensional 
procedure of section 1.4. and to either method 2 or 3 of 
section 1.2. 

As mentioned above, it has been shown from the trans- 
formation properties of the Hamiltonian that there are solutions 
of the Schroedinger equation consisting of mixtures of S and D 
states. These solutions are proper functions of the Hamiltonian, 
Hj the total angular momentum J; its z-component ; the total 

spin S; and the parity operator, P. In the scattering of 
neutrons by protons the incident plane wave includes both S 
and D waves in which these operators have the same proper 
values. The amplitudes, £s y of the outgoing S and D 

waves are related to the amplitude® , o*p of the incoming 

waves by a unitary, symmetrical scattering matrix S. Thus 
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(5.1) <• S » oL 

where 




We know and ^j). To find the differential cross section 
we must determine the elements which give ^ and ^ i) 

by (5.1). Actually it is sufficient to determine two of the 
( complex ) elements S,;'. since due to its symmetry and unitary 
character, the matrix S involves only three independent 
constants. 

Consider the expression 

(5.3) V V ’ (^) , 

l»\ ivl 

where and v^/ are the solutions of the Schroedinger 

equation, 

,5.4) 0, 

which have the asymptotic form 



39. See ref. (18). 
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4- t 
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^-*KA 
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4* 
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Here a»,a, describe the normalized angular and spin 
functions of the S and D wave respectively. Now let us use, 
in (5.3), trial functions whose asymptotic variation is 



O 

0 ^ 






(5.6) 



Then 



(5.7) 



S ^ 



in 



O'* e * >° l 
5 -A 

(») „*** 



ii s tjn. 



e 









Six, - 



l,) 2 - % t (v> _ 5 9 



Du 



Du 

(I ) r it* 









Hence, to the first order, we have 

i« / . 



(5.8) 



-r 2 vt f 11 ; a 11 «• > „ (*• n 

lx, -A -T- L *s -* <*t> J 

i« r "* >xi ' i) a u, i 
6 ^ 1 < X S ■+ <*0 1 ; 

or, in the matrix notation of (5.2) 

*r , U) 

4 _ ol s 



(6.9) 



Z K ,HI -X,) 

** ^’’ S V\’ 
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This enables us to obtain all matrix elements j . For 

<o u/ 

example, to find. Jjv we use trial functions with o^p =-0 , 

* V, -= O in which case the right hand side of (5.9) becomes 



K ( <•' x [ $»• \ / 0 \ U) {i J $ 

f I 015 °Hs„ sJUjV = T ef ‘ 0,0 



(5.10) 



and an application of the Rayleigh-Ritz method will yield the 
required matrix element. 

Let us note here that in this way we obtain a stationary 
result for the differential as well as for the total cross 
section ( The stationary character of the differential cross 
section is lost in the method of proper phases, described below.) 



The Structure of the Collision Determinant. 



Suppose wd want to calculate the element by the 

Rayleigh-Ritz method. We use as trial functions 




40* For calculating we noed not use different functional 
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to i* j 

where the radial functions U. and hr vanish rapidly as 

=? c>*=> . Thus each trial function consists of an S and 

a D part. Using the functions (5.11) in (5.9) we find that 

the tensor character of the potential introduces elements 

(•) , to 

connecting all coefficients of f* with those of 

The unknown is introduced by the right hand side of (5.9) 

of m 

and occurs only as factor of • 

On using the stationary property of %\\ , we are led 

to a system of equations precisely as in section 1.2. These 
are compatible only if the determinant of the coefficients, 

* vanishes. In the present case Ai\ has clearly 
the following, symmetrical form: 




, to 

forms for 'V* and Mr* . We use them here because for 
calculating off-diagonal elements of S, different functional 
forms obviously aro required. 
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, v I 

Because of its symmetry Aq > where is the 

large upper right hand block of Aq ♦ is therefore 

determined by the linear equation C.O . comprises 

two blocks of elements connecting S and D functions. These 
describe the interaction between the S and D state. In the 
limiting case of vanishing tensor forces they become zero 
and 



(5.13) 



assumes 


the 


degenerate 

s 


form 

D 






f" * 


0 


0 


A 1 


'X. 


% 


D 


0 


— 5? 












0 


0 


* 






0 


0 


X 


7t 



Hence is determined from the equation 



(5.14) 




whose elements no longer involve any D functions. This 
expresses the fact that the S wave is now decoupled from the 
D wave. In fact (5.14) has the same structure as (2.14), 
which was arrived at by considering only the S waves. 
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(r) 



Proper Values of the Scattering Matrix . 

The unitary matrix S has two eigenstates and 

with the property that 



4 



*■) 



(5.15) £ £ ‘c/ A) ~ 



A«) 



z- x 7 r. 



The X are the proper values of S and, because of its unitary 
character, may be written as 



(5.i6) 






PL _ A* 



Thus the eigenstates have the asymptotic form 



(5.17) 






.ID 






(1) ( 
Oi s 




Z ^X 

- e 1 




* V 


A 






«p w ( 


A 


^ Is. 
- e 1 


•A 



- - 1 



a e 



1* 



)n- s 

. s' w 4 'jj-'j (V s ~ 5 + D-d) 



(5.18) 



and 



^ ^ * T 1 * . Si~(Kflt (dj^iLs + c/^IIp) 



are thus standing waves with real radial functions 
and for that reason they are simpler to deal with than other 
states which are, in general, complex linear combinations of 
<f and Cj^"^ . Furthermore Schwinger ^ has shown that the 

41. J. Schwinger, unpublished lectures. 
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total cross section of the S + D mixture of a plane wave can 
be expressed in terms of the two constants Vj^. and only, 

namely 

(5.19) 6 ~ ) 

It i3 therefore of interest that tan!|^. and tan can be 
directly obtained by our method. In fact this problem is very 
much akin to that of finding tan in the case of central 
forces ( section 1.2. ). 

Y/e could proceed from the most general equation (5.9) 
but the following approach is slightly simpler; 

The correct 



II) 

C|> and 






both have the form 



(5.20) 



v{ 



<X, ^ UK, 






A 



A 



-V | 



<1* 



■+ 



l. 



UJ 



* 

y-) 






) 



where the a . , b ; , u* and w* may be taken as real and 

A 9 * A ' a * a v T x ’ 

Consider now the expression 



(5.2D x ^ y. 1 - v) ^ (<* r) 

where we admit only such trial , whose variations have 
the asymptotic form 

Sl i + 












D . 



(5.22) 
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Proceeding in the usual way, one finds 




(5.23) 

Therefore, corresponding to (2.23), we find 



(5.24) r + K (<h <*), -+ 1, Xz, ) ^ k ( a ^ Xou. 



as a stationary expression for tan . 

When we substitute a trial function of the form (5.20) 
in (5.24) we obtain an equation of the type 

(5 ' 25) ^ ^ ** «•') + *„• 

Proceeding in the usual way to a set of linear equations and 
their compatibility condition we find the following 



(5.26) 



Cl - 



for tan h : 
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Clearly this is a quadratic equation for tan ^ , yielding 
approximations for tan^ an< ^ ^ an * Th® ®l arnen ts 

arise from the interaction of S and D functions. In the absence 
of tensor forces they vanish so that ^ factors and we are 
led back to two independent problems for the S and D wave 
( compare (2.14) ). 

Prom the two solutions for tan we can find the total 
cross section by (5.18), whose error will be of the second 
order if our trial w ave function has a first order error. We 
can also find the differential cross section to the same order 
of accuracy as our trial function by solving the linear equations 
derived from (5.25) for the d,; and ^ * , thus obtaining the 



wave functions 






(i' 



and 



a 



and then taking that linear 
combination whose incident wave agrees Y/ith the incident S and 
D waves contained in a plane wave. To obtain the differential 
cross section to the second order of accuracy, we must fall 
back on the more general but more complicated equation (5.9). 



numerical Illustration . 

By using the method of proper values of the scattering 
matrix, we have calculated tan ‘IT and tan '■jf for neutron- 
proton ^ Oj scattering at 2.82 MEV. We have used 

the constants of Rarita and Schwinger * to have- a comparison 

Ml 

with Schwinger's numerical integration for the same energy. 
Thus the following interaction potential was used: 

(5.27) V = j(r) + 1 ’ lx $( r ) j 
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where 



(5.28) 



1- 




f(') " 


- t ; * \ 




0 } A >» 


II 






O ^ ^ ( 



and 9,^ is the spin dependent tensor operator. 



(5.29) 



< _ (fr * )C*i-r) 



6-,.^ 



A 1 



-13 



If is measured in units of =- 2.80 x 10 cm. 

Schwinger^ has shown that the Schroedinger equation 
leads, in the present case, to two coupled equations for the 
radial functions U and W r (- / ') of the S and D wave, viz. 



(5 






(5.31) 






u nN . 



where ^ and ^ are defined in (5.28) and k- 1-* 2. ' ^ 
is simpler to base our variational procedure on these 
equations, rather than to go back to the original three 



It 
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dimensional Schroedinger equation. We define 



( 5 . 32 ) 



U) 






U 4k k) 




if we use trial functions with the asymptotic form 



( 5 . 33 ) 



we find 



( 5 . 34 ) 



^ > $ S ( ^ H: “tov. Vj/ ; Oen IO\) 

^ ^ ^ D ( Ka ~Y V|^ c^-7 \4 a| 

V. ^ gU J \ 6(a 

- — K. ^ A s 5yt<X^ Vj ^ £ AaL^ ^ J 









Now the exact solutions which we are looking for have tan 
- tan z. tan vj . Hence the equation 



( 5 . 35 ) 



X-t J 

- K ^ A s 4 ^d) Ttcu^ vj ^ 



wherd n t" indicates the trial function as usual, gives a 
stationary expression for tan ij • 

We have used as trial functions 
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W 

H n - X e 4 A 4 

(5. 36) 1 

*X ^ 

I 

Id. 

and tan ^ were determined from the 
requirement that at A ^ 1 the functions (5.36) pass smoothly 
into the free particle functions 




K- lj 2-, 3 



Jl< ) 



^- 1 ^ 2 . 



(5.37) 



~ ^ Ka V|^ , L^yo \CA ^ 

* =-- A i>(. _ hi 104 ) 4 -A~ v u ) wi) 



A > | 



As checks on our calculation we have at each stage computed, 
by the variational method, the S and D phase shifts when tensor 
forces are absent and compared the results with the exact 
values which are easily obtained for the uncoupled states. 

The results of our calculations are collected in Table II. 



N 
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Table II . 

S + D Scattering with Tensor Forces* 



Trial Functions 


Uncoupled 


Coupled 




r 

O 

-4? 






■x. 

U S C | A + c i A, J 
w - d , a' ) 


6.9850 


- 0.3199 


6.986 


- 0.0003 


'u 

u - c f A + c 'j 
w » d ^ A* + d r A** j 


6.9850 


- 0.3070 


- 1.525 


- 0.0002 


X. 1 

U C^A + C^A + 

W d ( Ar^ + dtA* 1 j 


9.1398 


- 0.3070 


- 1.425 


- 0.0002 


Correct Functions 


9*2314 


- 0.3069 


- 1.435 


0.0000 



It is seen that the last pair of trial functions gives 
good agreement with Schwinger's numerical integration. 



1.6. Neutron Deuteron Scattering . 

So far we have been concerned only with two particle 
scattering processes*. We shall 3how that the methods which 
we have developed can be extended to many particle collisions* 

In this section the simplest collision of this type, namely that 
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between a neutron and deuteron, will be considered. Our 
procedure will be modeled after method 3 in section I. 2., 
but method 2 is equally well suited for many particle 
collisions . 

The wave equation for the system of the two neutrons, 
"l" and "2", and the proton, "3", is by (0.2) and (0.5) 

(6.1) ^ T «• \/( l i) + (/[»?) -t vj 

where T can be written as 



(0.5) ;T r. 




To remove the cumbersome factor in T, we make the following 
definitions : 



\ 




jr \ 




4M 


v( ; i) 


d 




E " 


V'pJ 




, £d / 



• 42. For the definition of the relative coordinates, see (0.3) 
and (0.4). 
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Thus 



(6. 3 ) D -- - g f ;) 



We now define 



(6 



.4) ?- ^ -ir ^ 



A £ 



t> 



so that (6.1) can be written in the form 



(6.5) 



T-'JiO 



Now consider the expression 



(6.6) I a => \ 



Y, 



) 



where V ^ is the volume bounded by the five dimensional 
surfaces. 



C ( " 

2.0, : * 



(6.7) 



l 



(v) 

A : ^ 3 



Clearly 

(6.8) 



l/V 



r ,J (f* ^1 ‘ 



Y/e aim to convert this into a surface integral over the surfaces 
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2. and. 2-<x and therefore express D symmetrically in 

terms of |r, and , namely 

(6.9) P - - ( V- V Ki . f ri + Vr^) 

as can be easily checked. Hence sr is the sum of three 
terms. 



(6.10) Si ^ R -f L+h 

where 




$ V ^ (A)^) 

Va 



(6.11) L 



[ (** V'* t* 4fV v 7 ri f) Wk , 

Va 



m - ( - K)(^) 

According to the exclusion principle Mr must be 
antisymmetric al in the indices 1 and 2. Hence if we insist that 
our trial functions have the same property, then will also 

be antisymme trical so that M - K. K and L can now be expressed 
as surface integrals by Green’s theorem: 
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(6.12) 



k = ( \ (. 'i y f - si w, 'i *) 

d ( (* 



o w <« 



X Jf- Z*.V Xx <Sr*)? 






where T, - *’| j A j , O' ^ is the three dimensional space ^ 4. 

and 5 ^ is the two dimensional surface .A, - Q. . Thus, 
the integral (6.12) extends over a portion of the five 
dimensional surface '2.^ * Similarly 

L = J ^J^V ri S^-SfV ri $ (*r v ) 



(6.13) 






, ( U* ** ) ? *, (<KJ 

# q> 

since, because of the antisymmetry of 1 , the two parts 

of L are equal. 

We now insist that a is so large that the follovying 
two conditions are satisfied: At /tj ^ , the wave function 

has split into a function ^(.^i ^describing the deuteron, 2 and 
3, and a function %■ (^describing its motion relative to the 
neutron I • i.e. 



(6.14) 



^ -- ‘{UO F(, t ) 



a ^ a 
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where ^ represents a negligible contribution from decaying 
modes. And further a must be so large that 



(5.15) 



negligible. 









2 ^ 

3 



^ V 



We restrict our trial functions to those which satisfy 
corresponding conditions, viz. 



(6.14') 







and 



(6.15») (ft) => negligible, f, >, | 

We assume that both «€(fJ and *?*(«) are normalized. 
| will have different forms depending on which phase we 
are dealing with. We shall assume that its angular part is 
normalized to Mr tr . Thus for S - scattering. 



(6.16) 



*r(r,) = ft ^ 

' -b 

6 (ft = -%a~ >| 



and similarly for higher angular momenta. 

If we neglect ^ and , the integrands occurring 

on the right hand sides of (6.12) and (6.13) are expressed 
most simply as functions of ft and ft ( see (6.14) and 



I 
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(6.14') ) so that it is of advantage to use these coordinates 
as variables of integration. Therefore we transform from 
ft . ft to ft • fv according to the equations 



(6.17) 








which are a consequence of the definitions (0.5) and (0.4). 
It follows that 



(6.18) [etr, ) (otft ) =• ) 

Also:' 

^■)ft %)r x * 

(6.19) 

Now for /i ( - a » ^ has the form (6.14) and 13 

(6.20) Sf-c^p,) F(r,) . 

Hence the gradient terms In (6.12) are 

( 6 . 21 ) | cf'(?i) f(ft)?, -» <{(?,) PV.) ft 
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V* -- I s ^ r ') i-* 5 f '(*,)*, 



( 6 . 22 ) 



+ ' F^») f + 5 <{(?.) F '( r .) ?, , 



where % r, i f, /if, 1 and p'/Vi) ) i F ^ r 'j 

denote the derivative of Fk) and along the outward 

normal. Let us first integrate over the ^ y -space, keeping 
Tj fixed. The integration extends over all points with the 
fixed and Ft ^ . Since c. Q. , it follows from 

(6.17) that at all points outside this region ^ ^ ~ “Q. ; 

hence, in virtue of (6.15) and (6.15' ), we may extend the 
integral over the entire -space. Furthermore, since <{(?<) 
and S£ (l?i) are independent of the direction of we can 
first carry out the integration over the angle between IT, 

A A 

and <? , in which all terms involving ft * give zero, 

A A 

while those involving give Mu • Hence 

cO 



^ * (30 ^ **- s i ( 4 " 

ii\=° 



< 1 ‘ 



(6.23) 



F y (0 JF'M 

-/(fO P v '( f i) ^ F ( r i) -<{l?i) 5 <(to) F ^') F Yn) 



\ ^ i ^ F ' 

V iio 
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Notice that 






due to the normalization of 
does not appear on the surface. 

Let us now specialize to the case of S - scattering 
and assume that for -*,?• a. , f* (<■.) has the form 



(6.24) ^ 






-V 



6 



) 



v/here is chosen to give a smooth logarithmic derivative 

across has the form (6.16). Hence 

<jf( a) =. <J b . 

(6.25) t\ 



jF'(a) = -<3B K 



'“a - 



so that 



(6.26) 



R 



- ■ Mu K. 5 ft 'A • 



Turning now to L we observe, by (6.19), that V f acts 
only on functions of l?«l and that 

(6.27) \j Xt . [ fiX' *{ I?' 1 ) ?>• 

A ^ 

Hence the integrand contains P. * as factor and the 

\ A 

integral over the angle between and vanishes. 

Therefore, by (6.10) and (6.26) 
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(6.28) 51^ - 2- ^ 



TR A S £ 

<? 



so that 




Zl 



u k. -Aft 



A 




(6.29) 




is a stationary expression for tan ^ • 

In the derivation of (6.29) v/e have met some difficulties 
which did not exist in the two particle problem: The neglect of 



^ and ^ in (6.14) and (6.14* ) and the assumptions (6.15) 
and (6. 15') introduce inaccuracies for any finite value of a, 
even if the interaction potential ha3 a sharp cut-off ( square 
well ); this is in contrast to two particle collision, where 
the wave function has accurately the asymptotic form for all 
distances larger than the cut-off distance of the potential, 
and where assumptions of the type (6.15) , (6.15’) do not 



This difficulty can be avoided by setting q. - „o but 
then another difficulty is encountered. We have in our 
treatment neglected all second order terms. However, unless 



(6.30) f*(V t ) 0*. f K 



enter. 




has the asymptotic form 
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where <fW 1 s the exact deuteron function and is 

exactly sinusoidal with wave number K- » the second order 
terms are appreciable over an infinite portion of the 
configuration space, and their integral, far from being 
negligible, makes infinite. Thus suppose that our trial 

function tends to 



(6 



•si) | f ($>.) l ^ r i 



instead of (6.30)^ where f* is sinusoidal. For large , 
(6.32) p -5> (y ^ jf M** ^ -y " £ p j 



where H y is the deuteron Hamiltonian, 



(6.33) ~ 4 

Hence, from the region Jlj > , say, in which (6.31) and 

(6.32) hold very nearly, gets the contribution 



( 0,0 

ax - ( 

\ A 




) + Sc^)] *jr*[r x ) 



(6.34) 



(-W &4 £ ^ ( r O^O(*iJ 









Of 1 f (a) i v ( <lr l J • 
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where is the expectation value of H -p in the trial 

state Cfy + 3 c^p . Clearly &X although a second order term 
( since is stationary ) is infinite. 

We are thus faced with the following dilemma: 

(1) If we use only an approximation for we must 
take a rather small a or else the effect (6.34) which is 
infinite for infinite a becomes important and our neglect 

of second order terms is unjustified. However, the finite a 
introduces an inaccuracy. 

(2) If we take a infinite we must use the exact deuteron 

function in the asymptotic form of which may be inconvenient 

for practical purposes. ( Such a procedure would bd nearly 
impossible for more complicated collisions. ) 

The way out of thi3 dilemma is suggested by (6.34), which 
shows just what the infinite contributions of the second order 
terms are and hence permits us to subtract them. This is done 
by replacing X<\ , (6.6) ; X^j Xq*# defined by 

(6.35) X cl - ^ ^ 

where 

?'* 



? 4 fr o - <£■ p 

Vf, •< Y | “I’ll - Uu -D)| -f K 4- 6 o . 



(6.36) 
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Thus I<\ depends in two way3 on the trial function: 

(1) Through P whose form is determined by the asymptotic 
deuteron function. 

(2) Through the trial - ^ itself, when substituted in 
(6.35). 

We shall show that 



(6.37) 1 j, v 12 ir k AB* - -nK AS ■+ iS ^ 

^ 2 , *3 ) 

. V 

where Q is the integral of second order terms which decay 
exponentially in all but a finite region of the configuration 



space. 

Let us denote 



\ (oh, )(Ar t ) 

V a 



(6.38) 


i P'lf 3 


Then, by 


(6.36), 




1*1%) - !*($*) 


(6.39) 






-* ) ( 6 >' - 


Now 


V* 

* V 


(6.40) 


X (\ s o 



3V(<fJ 

+ (c 0 ' - C») l | ^1 (olr, ) 

Vo. 

$J*(f) +-ld s *) 

&) 
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and, as wag shown above, without neglecting second order 
terms we have 

(6.4D 

where lim c 0 
« •*><£> 

To estimate the last two terms of (6.39), we decompose 
into 

v/ 0) 

V * } < & < °- 

(6.42) V a : ^ L > M 

n K) * L < a ) 

(X) 

where b is so large that in V ^ ^ and ^ ^ have their 

asymptotic forms 

(6.43) $ -<{(?,) F(r,) , ^ - c^(?,) + *K) 

or the corresponding expressions with 1 replaced by 2. 

Clearly the last two terms in (6.39) are both of the second 
order, so that their integral over the limited is of the 

type ^ . 6n the other hand in we shall find that the 

last two terms, give zero owing to the integration over fi 
( or o x ) . Thus in the region of V ^ where H, >Jb - , they give 
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-♦ (<V-e D ) 

+ (5=0-63.) Cf/l ^\{d^) 

= ' ^ r 3 ‘f* (“ W P + &P ) (f-t (« ( f, ) 

= 0 

Hence 

(6.45) [5 'i ) + - <C 0 ) (Sf^ j r: & ^ 

Prom (6.39), (6.40), (6.41), (6.45) we now find 

which is the statement (6.37). 

This means thati&n^as calculated from 

(6.46) +|ltK A fc* 5 U (j 

is stationary if we admit only such trial functions whose 
asymptotic form is 

(.<?,) ^ ft -t 6* 



(6.47) 
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Thus the exact form of the Internal wave function and the 
binding energy of the colliding nuclei are not required. 

Numerical Calculations . 

We have made exploratory calculations on n-d scattering 
in preparation for a more extensive investigation. For 
simplicity we have so far only used trial functions of group 
structure, but vie remark again that this limitation is by no 
means inherent in our method. 

'We have assumed an ordinary Wigner potential 



These constants give the correct binding energy of 2.18 MEV. 
for the deuteron. 

As trial functions for S-scatterlng we have used 




with V 0 — - 72.00 me’* 



HI.) 

and a 0.8 in unit3 of A 0 =. 2.80 x 10 cm. 



(6.49) 





where vie have set 





43. These constants are almost identical with those of Motz and 
Schwinger, ref. 23, who however used an exchange potential. 
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(6*50 cont'd) 

4 ) db) dli) 

1U = ~ f > ii')Ci 

All distances are measured in units of 4 q . CJ^ ( ^>J is the 
best approximation of its type for the deuteron function. 

The polynomial approximation for {'(r) has the wrong 
asymptotic form. Hence it was assumed to hold only inside the 
region A < <\ . 

Y/e have calculated only S-scattering at zero energy. Hence, 
asymptotically 




where A and B are chosen so that the polynomial form of 

passes smoothly into (6.50*) at As ^ . 

As we are limiting ourselves to a finite region, the 

difference between 1(a) and I* (a) is negligible ( see (6.35) 

and (6.36) ). The calculations to be described were carried 

out by using I* (a) throughout. 

Y/e have calculated limit ( tan \/ ) =. X, for quartet 

X •* o 

and doublet scattering, which determine the total cross section 
in the limit of zero energy by the equation 



Sc-oJ$3A/ 






PouJLfot 
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(6.53) <f 



Hu (J ^ 



where q and d refer to quartet and doublet scattering 
respectively. 

In the case of doublet scattering, the convergence 
with increasing values of V\ was slow. Since in the case of 
two particle scattering method 2 gave somewhat better results 
than method 3, ( see Table I, p. 30 ) we have calculated X ^ 
not only by (6.29) but also by the equation 



which bears the same relation to (6.29) as does (2.32) to 
(2.23). X is determined from the logarithmic derivative L 
according to the equation 



However, the convergence was only slightly improved. 

The approximate values of X, obtained by our variational 
methods depend both on a and n. For large enough n we expect 
the following variation with a: 

(Al) When a is too small the expression (6.29) for tan 
is incorrect. Hence X should depend strongly on a. 

(A2) When a is sufficiently large for (6.29) to hold, 

X should be fairly independent of a until 




8 



(lf%) - ?'(*)) = o 



(6.52) 
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(A3) a becomes so large that F n can no longer 
approximate the true function over the entire interval. Our 
variational principle than breaks down because second order 
terms become important. Here again we expect a considerable 
variation of X with a. 

Thus we expect to find an extremum in the region (A2) 
and the value of X at that extremum may be taken as our 
approximation for X. 

As we increase n we expect 

(Bl) Gradual extension of the "plateau" (A2) towards 
larger values of a. 

(B2) Increasingly smaller differences X ( a; n + 1 ) 

- X ( a; n ) for large enough, fixed, a* 

In an approximation method n should be increased until 
(Bl) and (B2) are so well satisfied that further increase of 
n is not likely to change X appreciably. 

In virtue of the choice of our potential and trial 

function, all integrals occurring In 1^ could be analytically 

are 

evaluated. Our results for X<^ and X^ A given In Tables III, 
IVa and IVb and plotted in Graph 2. The best 3 - parameter 
wave functions ^ [redetermined in the usual way from our 
variational principles are shov/n in Graph 3. 



GRAPH 2. NEUTRON 8EUTER0H SCATTERING 
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Table III. 



Quartet Scattering : 




(ajn) 



/^yp » 



( Method 3, (6.29) ) 



0 , /AoS^ 


0 


i 


2 


2.5 . 


- 4.39 


- 2.16 


- 1.58 


3.0 


- 4.39 


- 1.91 


- 1.85 


4.0 


- 4.39 


- 1.84 


- 1.83 


6.0 


- 4.39 


- 2.03 


- 1.84 


8.0 


- 4.39 


- 2.26 


- 1.92 




Table 


IV a. 




Double t 


Scattering : (ajn) . 




( Me thod 3 


, (6.29) ) 






0 


1 


2 


2.0 


+ 8.53 


- 1.18 


- 1.58 


2.5 


+ 8.53 


- 2.03 


- 3.73 


2.6 


+ 8.53 


- 2.14 


- 1.03 


3.0 


+ 8.53 


- 2.51 


- 1.58 


4.0 


+ 8.53 


- 3.38 


- 2.01 


6.0 


•h 8 * 53 


- 5.65 


- 2.98 


8.0 


So 53 


- 9.16 


- 4.05 
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Table IV b. 



Doublet Scattering : (a;n) / A 0 . 

( Method 2, (6.51), (6.52) ) 





0 


1 


2 


2.5 


- 3.53 


- 0.36 


- 1.51 


3.0 


- 4.62 


- 2.04 


- 1.47 


4.0 


- 7.53 


- 2.85 


- 1.90 


6.0 


- 20.23 


- 4.47 


- 2.78 



Interpretation of the numerical results: 

A. The convergence of X^ is much more rapid than that 
of X ^ . This is due to the fact that F( f ) is more irregular 
for doublet scattering ( see Graph 3 ) and therefore requires 
more terms for adequate approximation. 

B. The variation of X ^ with a and n is as expected 
( pp. 73 and 74 ) ; X ol seems to behave in a similar manner 
although the number of parameters is not large enough to 
exhibit a proper plateau. 

C. To determine the total cross section we must select 
the best values of X and X d from our tables. In the case 
of X ^ ( Table III ) we have evidently come very close to the 
limit which can be reached with group structure form. We see 
that our results depend markedly on a for values smaller than 3* 
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This indicates that the interaction is appreciable out to 

about that distance* We therefore choose X ^ =. X ^ (4;2) - -1.83^\ 0 

The choice of X ^ is more difficult. At a = 3 we cannot be 

sure that the interaction is completely negligible; on the 

other hand at a r. 4 we are still, quite far from convergence. 

Let us nevertheless choose X ^ - X ^(4;2) - - 1.90 A 0 ( see 
Table IV b ). This gives, by (6.53), 6 - 3.39 barns. If 

we estimate the limit of X ^ which could be reached by taking 
larger n, as X ^ •=. - 1. 50/i^we obtain ^ =- 2.94 barns. These 

figures are in qualitative agreement with the experimental 
value of 3.2 barns 

Let us summarize the possible inadequacies of this 
tentative calculation. 

(a) Wigner forces may not be adequate to describe this 



process. 

(b) Our approximation for ( see (6.50) ) may 

not be sufficiently accurate. 

(c) The group structure approximation may not be 
adequate. 

(d) The Schroe dinger theory may not be capable of an 
accurate description of this process. 

Of these shortcomings (a), (b) and (c) can be remedied 
by more extensive calculation of the same kind. When this is 
done, it will be possible to decide to what extent (d) Is true. 

In future calculations the following points should be 
given consideration: 

(1) Terms allowing for the polarization of the deuteron 
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should be introduced. 

(2) It may be of advantage to use trial - F's of the 
correct asymptotic form. The reason for this is the following: 
If we use a polynomial expression for F (*'’), we must include 
a large number of terms since, due to the loose structure of 
the deuteron, the region of interaction is very large 
(cs 8 x 10 cm ) . This necessitates the calculation of a 
great many matrix elements, which is probably more laborious 
than the numerical integrations which become necessary when 
functions of the correct asymptotic form are used. 



1.7. The General Case of Composite Collisions . 

The method used in the preceeding section can be 
extended to all collisions of composite nuclei in which triple 
or multiple disintegrations are energetically forbidden. The 
analytic aspects of the problem are quite analogous to those 
encountered in n-d scattering. However the combinatorial 
side of the general problem is somewhat complicated and no 
detailed description of the general case will be given here. 

We shall merely mention a few points of interest. 

Modes of Collision . 

Our picture is the following: V/e have a number of pairs 

of nuclei colliding and a number of pairs ( possibly different 
one3 ) emergin * after the collision. The groupings of nucleons 
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into pairs are alternative (not independent) processes, 
each pair comprising the same nucleons. There is, then a 
"region of interaction" in the configuration space, where all 
nucleons are close together, and an "asymptotic region", 
where the nucleons have separated into pairs. 

On a hyper- surf ace, , separating the two regions, 

the wave function ha3 the form 



where, if and describe the normalized internal wave 
functions of two colliding nuclei and F their relative motion. 



(7.1) 




ii) 



(7.2) (Q . 




») V>l 



A function of the type (7.2) will be called a mode of 



collision. 



is a proper function of a complete set of operators 
( H, J, P, A^..), its modes satisfy the following 




asymptotic orthogonality conditions: 





44. The term "collide" will be used in the remainder of this 
section to describe both collisions and disintegrations* 
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For if 



and $j 



describe different pairs of nuclei. 



they exist on different portions of the hypersurface 



Z 



and 



hence (7,3) is satisfied. If however they describe the same 

ii) 0> U> l l > 

pair of nuclei, then either or 4 , so that the 

orthogonality (7.3) follows from the integration over the 



internal coordinates of the nuclei. 



Number of Arbitrary Constants of the Scattering Matrix . 

Vie have remarked in section 1.5. that due to its unitary 

character and symmetry, the two dimensional - scattering matrix 

involves only three independent constants. In the present 

case the scattering matrix has the same order, n, as there are 

different modes of collision. When represented in the scheme 

of incoming and outgoing modes, we have proved that the number 

of independent constants is (n) (n+1). Only this many constants 

12 

need to be determined. 



Raylelgh-Rltz Method . 

We have already applied the Rayleigh-Ri tz method to two 
problems involving two modes of collision: n-p scattering with 
tensor forces ( section 1.5. ) and n-d scattering ( section 1.6. ). 
In section 1.5. we have discussed two different variational 
methods: One dealing with complex ( running ) wave functions 
and yielding all elements of the scattering matrix; the other 
dealing with standing waves and giving only the proper values 
of the scattering matrix. ( Only this second method was used 
in n-d scattering. ) 
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These methods can be generalized as follows: 

Running waves: 

When dealing with a special case of running waves in 
section 1.5. , we introduced a wave function, sr * > 

(5.5), whose angular and spin functions were the complex 
conjugates of those of the function 'tv This was done 
in order to make use of the orthogonality relation, 
corresponding to (7.3) 

In the present more general case, such a procedure 
is not convenient, since the wave function does not everywhere 
separate into radial and angular parts. On the other hand, 
if vie take the complex conjugate of the entire function, the 
sense of the incoming and outgoing waves is reversed and the 
oi ’ s and ' s, denoting their amplitudes are mixed up. To 

avoid this difficulty we use the following notation: 

(7.4) * VJ4/0 * 



where ^ 

> 






and 



^r> 



have the asymptotic forms 



* 



(7.5) 



> sr M f 

- 1 ^ p (.'h 





(A ) £ A 


* [ 




Jr*) e/"^ 

s* — JL 1 




J\. + l 


r -% 



H 



S * Is the angular function of relative motion of the mode 
We note that by the definition (7*4), 

< 7 . 6 , 



. Du 



$• . 
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r l*'/*) 

which is the reason for introducing i , as only in this way 
can the o L f s and ^ 's remain associated with incoming 
and outgoing waves respectively. 

We now form the expression 



.7 ) = $ ^* J ** 

V 



(7 



0> l,; 

where tL and J- have the same quantum numbers, V is the 
region of interaction, H the total Hamiltonian and E the 
total energy of the system. The first variation of 
can be transformed into a surface integral in the usual 
v/ ay and we are led to 



Si 




(7.8) 



tl ^-'>2 

I 3w ^ 

= Z 'i 

“ ' V " \ M; ; 



(Compare (5.7)) Therefore the following equation is 
stationary relative to variations of and 



(7.9) 




-'l 



_ _ C~') <■*> 

\L • d • 0 > . . , 





We would like to convert this into a stationary 
expression for the elements of the scattering matrix 
(see (5.9)). Now in the case v/here groups may collide 
with different kinetic energies, we must modify the 
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definition of the scattering matrix given in (5*1). The 
amplitudes are linear combinations of the OC £ , viz. 



(V.lu) 




The total flux into the region of the configuration space 
enclosed by £ is 



0- i k; \ K-r'i^r} 

(7.ii) -.1 ( S’; X- \ c - j l \otj 

} J 

•4iiK« ; r -zif| 

The last line shows tiiat the matrix 



(7.12) 




is unitary, and hence may properly be taken to be the 
scattering matrix. If in analogy with (5.2) we define 




46. A. Wintner, Spektraltheorie der Unendlichen Matrizen, 
(Leipzig, 1929) p. 34. 
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then by ^7.1u) and ( 7 . 12), 



(7.14) 'fe =. S • f\ 



In this notation i7.9) can be written as 

-r £ M " 11 ^ c A ( ^ 

,7.15) ^-z\ ■+ ^ ^ B Zfa^jL X ™ ^ 



Hence by choosing our trial functions with only a single 
non-vanishing component of A, we can obtain a stationary 
expression for any one matrix element of S at a time. 



Standing Waves: 

The modes of the proper functions of the scattering 
matrix all have the same phase shift (see section 1.5. ,(S’.lS')‘H. ) 
We form 

(7.16) 1 * S \P*C-w +e) $ 

_ V 

where ^ Is a trial function which consists, asymptotically, 
of standing waves. The first variation of X now becomes 

(7.17) ST ^ ^ \ 

where is the amplitude of SU A; jj\' A in F , the 

function of relative motion. 



Hence 
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( 7 . 18 ) 



I '+ ft; 2- -£a~ vj 



The collision determinant, which arises on application 

-tk 

of the Rayleigh-Ritz method is of the Y\ degree corresponding 
to the sum of n squares multiplying tan v|^ . Its n solutions 
are approximations for the tangents of the n proper values of , 
the phase shift. 
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II. THE VARIATION ITER A TION METHOD IN COMPOSITE COLLISIONS . 

II. 8. Preliminary Considerations . 

4 \ 

Schwinger has recently developed a variation iteration 
method for two particle collisions. His treatment leads to a 
stationary expression for tan , involving a Green’s function 
and giving very good results, even with rather crude trial 
functions. In this section we shall first rederive his 
results in a slightly different way, to show its relation to 
the bilinear variational principles discussed in Part I; and 
then shall examine this method with a view of generalizing it 
to many particle problems. 

Consider first the one dimensional Schroedinger equation. 



Let the Green's function .G(x, x') satisfy the differential 



( 8 . 1 ) 







- o } 



where u satisfies the boundary conditions 




equation 



( 8 . 3 ) 
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and certain boundary conditions to be determined later. 
Multiplying (8.1) by G and (8.3) by u, and integrating over 
x from 0 to a we find 





To make the integrated terms vanish we now impose the same 
boundary conditions on G as ar 6 satisfied by u, namely 



Q x l ) - o 

(8.5) 

6 i * ) x ) ^ ^ ) C ^ X, -V . £Vc } CiA X -> O0 ^ 

where F (x*) is some function of x 1 . The solution of (8.3) 
and (8.5) is 



( 8 . 6 ) 



c [*, *') 




S i/*' ^ ^ 




S ov* K-Af. £ ^ Kx 



In virtue of (8.4) and (8.5) we now obtain, on letting CL-^oO, 

(8.7) * 'k j G •{(*') u[V) ^ • 

0 

This may be regarded a3 a proper value problem for^-, 
for a given cot^^and the discussion can be completed from 
this point of view: (8.7) is multiplied by f(x)u(x) and 
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integrated over x, giving 

oo 



^ [M ^ = ^ ^ ■Sx^.K^C. ^(*) 



(8.8) 



ot> 



-1^ {I*) *A W ^ O l*y) deceit 



whe re 



(8.9) (*>*') - -1 Si’~v KZ- 



> 



(8.7) is a stationary expression for 7- corresponding to 
a given cot or for cot corresponding to a given . 

On the other hand we may derive (8.8) as follows: 



We note that u defined by 

L* 



(8 



CT-' 

.10) u L \ <?(*,*’) ^ 



J 



becomes , asymptotically 






( S i/w 14 



)• 



(8.11) 



X 4 ^(K 6^o |4* 

f si- nx' -f(v) uM JUc' 

K > 



and thus has the correct asymptotic form and phase shift 
for any trial u; furthermore by (8.7), u u when u is 
the correct solution. Hence if used in our bilinear 
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variational principle (2. 48), in conjunction with u ( s u, 
we have simply 



t> V 

(8a2) 7i( (£„(*, *0t ^ Si - -si- x*’) 

0 K 

° ') 

tlzx* 0 



\ 



JoC 



for arbitrary variations of u. This is of course 
equivalent to (8.8). 



Three Dimensional Formulation . 

Turning now to the solution of the equation 



(8.13) 







one can proceed in an entirely analogous manner. The solution 
of (8.13) corresponding to a wave of unit amplitude incident 
along is written in the form 



(8.14) 



"hu - 



A 






A k lr-r‘1 
4tt Ir-r*) 






and becomes asymptotically 



(8.15) LKj4 ^ 



^ e 



A U 2 .f 



; ra 









R -f 



-A uf rr 

Henco the scattering amplitude in the direction 1^ is 
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given by 

(8.16) ^ K,k.) | e' AS ' r {(>') ^ ^( r 'J ^r'J . 

This permits us to write (8.14) in the homogeneous form 



(8 - 17) W -- ^ jc(r,,>) {[*■; ^(r'j ( dr ‘) 



where 



(8.18) 







For large values of r, (8.17) becomes, in the direction ^ 

'('HilO ? 

4 

(8.19) 







-f l v ‘) Hr l ’■’ ) (« lr ' J 



Hence no matter what function 4^(y) 13 subs ti tuted 
in the right hand side of (8.17), the left hand side always 
has the correct scattering amplitude f( ^-i.jH.|). 

This, in fact, is the. reason for replacing (8.14) by (8.17). 

Returning to the variational principle, it is clear 
that if we use only such ^^which have the form 



<(*') . 



( 8 . 20 ) 
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I. / / 

and arbitrary then 3ince ^ (Kj,, K,) ^ '{(* 2 ., K,J 

It reduces to ;I „ -o* When this i3 written out in 

*“i K| " 

full, it serves as a variational principle for the 
determination of 



J. 

4 ir {(^.k,) 



t-K,(0 {[*) e'" tr ^dv') 







(8.21) -{(<■) Cj'njr) [tflv) 




(^,(0 f (0 



K (r-r'J 

I 



‘‘‘J yX ) (dr) (tl/'J -o 



Let us note specially that no boundary conditions had to 
be imposed on either vj/ ^ or ♦ This is in keeping 

with the fact that all integrals occuring in (8.21) extend 
only over the region of interaction, where f(r) and f(r’) 
are not zero. 



The one and three dimensional treatments have in 
common certain features which may be expected to carry 
over into a similar treatment of composite scattering 
processes : 

A* Both treatments involve one iteration^(8.10) 
and (8. 20 and hence give excellent results even with 
poor trial functions. In principle the accuracy can be 
increased by iterating a number of time3. 
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B. The trial functions need not satisfy any boundary 
conditions; instead, 

C. The boundary conditions are automatically satisfied 
by a proper choice of the Green's function. 

D. The Green's function 



the same differential equation as does the wave function 
asymptotically ( see (8.3) and (8.18) ) and 



amplitude ) as the correct wave function. 

II. 9. Neutron Deuteron Scattering: Six Dimensional Green's 

Functions . 

As a prototype of composite scattering processes we 
shall consider in this section the collision of a neutron 
with a deuteron. We shall generalize the method described 
in the last section in a straightforward manner, but shall 
find that for composite scattering it is too complicated to 
be of practical use. 

The wave equation for our system is, by (6.3) - (6.5) 



1. except for the inhomogeneous term, satisfies 



2. has the same phase shift ( or scattering 



(9.1) 




where 
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(9.2) 



D * -(yr^ + i \}£) - - 



M 

3 




Let us for the moment ignore the spin dependence and 
simply insist that the wave function be anti symmetrical 
in 1 and 2. Then, asymptotically the S-wave ha3 the form: 



r, : sir ( f(?0 1 i 



(9.3) 



The wave equation (9.1) becomes asymptotically: 
5 or* *. Lj - 



(9.4) 






: V' U,5 + € »)]^ 



According to D.I., section II. 8, the Green’s function 
should satisfy the differential equation 

(9.5) L £ 



in the asymptotic region of the configuration space, 
where L reduces to L and L for large values of n. and 

\ U. 1 

respectively. In the one dimensional problem L was simply 
3- + and made no reference to the interaction 

potential; hence a solution for G could be easily found. 
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The fundemental difficulty in the case of composite 
collisions is that parts of the interaction , namely the 
internal forces of the colliding nuclei, persist at infinity . 
Therefore they appear in the asymptotic operators L and L 
and hence also in L, which makes the determination of G 
practically impossible. In fact the simplest L which we can 
choose is L ; P + v lr> , since is the only interaction 

which disappears asymptotically. Of course the practical 
solution of 



is quite impossible. 

For this reason we shall not pursue this method any 
further, but merely state that if formally carried through, 
it has the properties A, B, C and D mentioned In the last 
section and leads to stationary expressions for the elements 



of the scattering matrix ( i and j refer to the incoming 
and outgoing neutrons, ^ and to their* direc tions ) . 

Another formulation which suggests itself is analogous 
to the treatment of the bound state: We write the wave 
equation as 



(9.6) D - 






•t y Q ■=■ — <£ ( (\~ tx) ^ 




A - 




and introduce a Green's function to satisfy 
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, 9 . 9 , (y r ; , | ) c(r 3;?3 j r,*, fi ') =-J (rj-fj'J J( fj - j3 'J 



A solution of thi3 equation is 



(9.10) 




where 



(9.11) 



■» = I U ' r; | , 



Since £p is negative in our case, G falls off exponentially. 

Hence, on multiplying (9.8) by G, (9.9) by * , subtracting 

and integrating over and ^ , the surface integral 

vanishes and we are left with 



(9.12) \F (q J £,') - ^ C (r^ + 

However the exponential vanishing of C- wipes out all reference 
to the boundary, from where in the past we had always derived 
the connection with the scattering properties. Thus the 
present formulation, while providing an iteration procedure 
for the collision wave function, does not lead to a stationary 
expression for the phase shift or scattering amplitude. 
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II. 10. Neutron Deuteron Scattering: Three Dimensional 

Green’s Function . 

1 Since, as we have seen, the straightforward generalization 
of the two particle procedure is of no practical use we have 
tried the following alternative approach: 

’Where f, , the wave function separates and therefore 
satisfies not only, ( see (9.4) ) hut also 

(10.1) r, -9 ^ '£($’,,1%) =0 

Hence, in partial agreement with D.2 ( section II. 1. ) we 
Introduce a three dimensional Green's function, which 
satisfies 

do.2) (y r ,\K v ) = -SK-f,') 



and has the correct scattering amplitude in the direction 
. Suppose that we consider the wave function, 
corresponding to a plane wave of neutron 1, of unit amplitude. 



ifk . 



incident along 1-4^ while both neutrons 1 and 2 are scattered * 
Then the scattering amplitude of neutron 1 along is 

* A Green's function satisfying (10.2) is 



46. This Is in conflict with the exclusion principle but 
convenient for our present considerations. By super- 
posing the symmetrical wave function, with a negative 
sign, a solution obeying the exclusion principle can, 
of course, be obtained. 
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(10.3) Hu In-/,'! 



• W Cl - 
a . r t 



i 

-x - r, 



It satisfies the differential equation and asymptotic 
conditions of the wave function only where 'C \ — ^ oo , but 
not where — 9> cO , and thus not over the entire asymptotic 

space ( see section II. 1., D.l. ). 

Proceeding nonetheless we write the wave equation (6.5) 
for our system in the form 



(10.4) (? 0 r.) 

where 

(10.5) - btC?o’ r .) H -iv Si l 4 U i3 4 »a - 

is a kind of effective potential, acting on particle 1 when 
is small. 

From (10.2), (10.3) and (10.4) we can deduce in the 
usual way that the wave function satisfies the integral 
equation 



(10.6) Sir K w [?„*,) * U(' r .)V-,’)Ulf,j»-;J l P^( ? , ) r l , )yr,') 



as may also bo checked directly. For large in the 
direction and small values of ^ this becomes 
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(10.7) 




Thus again, regardless of the ^ r/J used in the right 

hand side of (10.6), the left hand side has the correct 
scattering amplitude along . However unless 

i 

splits for large T, into a product of the form 



(10.8) ^ «\') — > <•{(?,) f (r,’J 



u*„ s > does not vanish so that (10.6) would diverge. 

For in contrast to is not a simple multiplicative 

factor and "cuts off" only wave functions of the type (10.8). 
Now consider the expression 



uo.9) x- j + 

Clearly, for the correct I r: 0. Let us now consider 
the variation of I, if we admit only trial functions of the 
form 



( 10 . 10 ) 



- 



G U r >') 5} uj 1 ' l fo r '.)(^ r 0 



where will be suitably restricted; and 

A 

unrestricted. Then 



I 3 

K <? 
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( 10 . 11 ) 



Si ~ [ ' k f (-0 }(°< r . 



As in section 1.6. , (6.8) - (6.13) this leads to surface 

integrals over the five dimensional surfaces ^ . * l A . - C\.) 

\ V N T ( 

and L ^ The one over Z. ^ vanishes because 

P / U) a) \ 

J owing to its expression in terms 



has the correct 



of a Green'3 function. The one over "vanishes if 
describes only outgoing waves of neutron 2, since x-K^ ; 
contains only such waves and hence the integrand vanishes. 

We shall show that has indeed this property 

( of outgoing neutron-2 waves ) if only has it. For 

A L ~ q and $ Q , exceeds 2/3 a. 

Hence, by our assumption. 



( 10 . 12 ) $ ^ f 




Now for ^ >/ 2/3 a, (10.12) is an exact solution of the wave 

equation (10.4) so that 

(10.13) [v/ 4 K % U(«f, Kj y = 0 

Hence * , 

' J G 1 r '^ r ') 

-- -\Q U, r,‘) 4 K 1 ) $kJ* C?«i *■') r >) 

- K) 



(10.14) 
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(10.14 cont’d) 



^ \ GU.r,') V r , [?„ 0 - 

$Ki" (f.,0 *r, 6 



But when 0 ^ S/3 a^and f, = ^ (. ft; ** J ” ^ / 



so that the surface Integral vanishes. Hence, using (10.2), 
we find 



t'UA, 



(io.i 5 ) r L 



<s>o 



•• ^k 5 * = # «-! ^ <?(?>) f e -x; 



as was to be shown. 

Therefore if the trial function d? satisfies the 

conditions 



^ ^ -=> ^ 




) ^ f l —5 c>° 



(10*17) 



47 K- r.'l 



^ <$nj'> (?„ r, ') (*tf, ) «•/ j 



4ir 



■{W’4 










(0 



; 
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provides a stationary expression for 

In this way we can obtain the amplitude for elastic 
scattering in all directions. In the simple case of n - d 
scattering, the amplitude for exchange scattering i3 then 
determined owing to the restrictions to which the scattering 
matrix is subject. Thus, let U3 consider a sub - matrix j 

corresponding to scattering with a given angular momentum, L 
where the subscripts refer to the incoming and outgoing 



J 



neutrons. There are two modes of disintegration so that, by 
< u 

section 1.7. ( ) involves at most three arbitrary constants. 

In fact, due to the equality of the two neutrons this number 
is further reduced to two. Thus a single diagonal element com- 
pletely defines the matrix. It is clear therefore that the • 
determination of 
be deduced is sufficient. 



from which all elements can 



Variational principles for can be obtained 

in a similar way, but in addition to volume integrals, as 
occur in (10.17) they also involve surface integrals. The 
same holds for variational expressions for the proper phases 
of the scattering matrix. 

No numerical application of (10.17) or its equivalent 
for the phase shifts has been made; The labour involved 
would be enormous. Thus the calculation for the phase 
shifts involves a six dimensional integral which can be 
reduced to five dimensions by an analytical integration. 

The remaining integral, owing to the presence of the Green's 
function must be evaluated numerically. Furthermore, since 
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the Iteration (10.10) extends only over three dimensions, it 
is not certain how much it improves the trial function and 
hence one cannot be sure that (10.17) will yield very good 
results . 



\ 
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CONCLUSION. 



The problem of many particle nuclear collisions is a 
difficult one. Even in classical mechanics the treatment 
of several interacting particles is in general complicated 
and hence one cannot reasonably expect a royal road to the 
solution of the corresponding nuclear problems. 

As remarked earlier, however, good results have been 
obtained for many particle bound states, with a reasonable 
expenditure of effort arid time. In Parts I and II we have 
examined the two methods which have been used for bound 
systems, for their applicability to collision processes. 

We believe that the Rayleigh-Ritz method, discussed 
in Part i, provides a practical procedure also for the 
solution of collision problems. To be sure, the fact that 
the trial functions must satisfy certain asymptotic conditions 
(e.g. (6.14) ) causes some complication. But the integrals 
to be evaluated in a collision calculation are of precisely 
the same kind as occur in the binding energy calculation of 
the corresponding bound state. This indicates that scattering 
cross sections can be obtained by an amount of numerical work 
comparable to that involved in bound state calculations. 

On the other hand,' the variation iteration method 
( Part II ), although yielding strikingly good results in 
two body collision problems, appears to be unsuited for the 
treatment of more complicated scattering processes. For •. 
we have seen that the theoretically most satisfactory procedure 



i 
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( section II. 2. ) is completely useless from a practical 

point of view, while other methods ( section II. 3. ) are 
still extremely complicated and of doubtful value. Of coursp 
one cannot rule out the possibility that some better iteration 
procedure exists. However any Integral formulation increases 
the dimensionality of the integrals to be found. Thus, unless 
some unforeseen simplification can be made, we do not believe 
that any such method can be as useful as the Rayleigh-Ritz 
method. 

While this paper falls far short from providing a 
complete solution of the problem of composite nuclear 
collisions, it is .hoped that it has made some contribution 
towards paving the way for a systematic treatment of such 
processes. 
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